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ABSTRACT 


A remote sensing technique is described which 
utilizes elastic scattering and rotational Raman scatter- . 
ing of laser light in the atmosphere to obtain soundings 
of turbidity, transmissivity and density. 

When irradiated by monochromatic light, both molecular 
and particulate components of the atmosphere scatter light 
at the incident frequency. In addition, the gaseous 
component produces a temperature-dependent pure rotational 
Raman spectrum which extends about 200 cm"* 1 above and below 
the incident frequency. It is shown that particulate scat- 
tering, as well as that from the minor atmospheric constit- , 
uents, is negligibly small in this interval, in comparison 
to that from N 2 and 0 2 . 

A scheme is devised whereby, through selective weight- 
ing of the rotational Raman lines, the effect of atmospher- 
• „ 

ic temperature structure may be eliminated. 

The close spectral proximity of the elastic and Raman- 
scattered signals, combined with the fact that the Raman 
scattering is quite weak, produces special requirements 
for the spectroscopic and light-gathering components of 
a rotational Raman laser radar system. These require- 
ments are investigated, with a view toward the design of 

xii 



a practical system. 

A computation of typical signal-to-noise ratios is 
made. It is shown that daytime signal-to-noise ratios 
greater than 10 db are to be expected for observation 
heights of 5 km and below. For nighttime work, 10 4b 
signal-to-noise ratios are achievable to altitudes as 

high as 15 km. 
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Chapter I . 

INTRODUCTION 

Several investigators have utilized optical probing to 
obtain remote density and turbidity profiles of the atmos- 
phere* Their studies depend on the fact that some- photons, 
in passing through the atmosphere, undergo elastic colli- 
sions with gas molecules and particulates, and are thereby 
scattered away from their original paths. Some of these 
photons are scattered toward the receiver and constitute 
the observed return signal. All of these studies have in 
common the problem of separating the return due to colli- 
sions with molecules from that due to collisions with 
particulates. Furthermore, since knowledge of the number 
of photons returned in a given direction in an atmosphere 
containing particulates does not uniquely determine the 
amount of signal attenuation or the amount of scattering, 
some assumptions must be made in this respect. As a result 
of the latter problem, most work in this area has focussed 
on stratospheric probing, where particulate effects are 
relatively small or confined to discrete layers, and where 
signal attenuation is nearly constant with height* 

Elterman (1966) assumed a molecular number density ^ 
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distribution and confined his study to determining profiles 
of attenuation due to particulates. He further assumed 
that particulate attenuation was zero at 35 km, except when 
that assumption led to negative particle attenuation values 
at lower altitudes, in which case the attenuation at 35 km 
was adjusted to the smallest value that would give values 
greater than zero at all lower altitudes. Since Elterman 
employed a searchlight system and determined his observa- 
tion heights through triangulation, his observations for 
different heights were made at varying scattering angles, 
and thus it was necessary for him to assume a phase func- 
tion, or directional distribution pattern, for the particle 
scattering. Had he used a laser radar System and monitored 
only backscattered radiation, it would not have been neces- 
sary for him to concern himself with the details of a phase 
function, but instead to assume some relationship between 
the number of photons scattered toward the receiver and the 
total number scattered out of the beam. 

Fiocco and Smullin (1963), Fiocco and Grams (1964), 
Clemesha, et.al. (1966), Collis and Ligda (1966), Kent, 
et.al. (1967), Schuster (1970) and others have' directed 
their analyses toward finding discrete aerosol layers in 
the stratosphere. In each of these studies, the observed 
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stratospheric profile is fitted to a model of the expected 
profile for a molecular atmosphere, and the small regions 
of deviation are ascribed to the presence of aerosol lay- 
ers. 

Kent, et.al, (1971), on the other hand, fit their, 
soundings to a standard atmosphere at the 30 km level, and 
assume that all deviations above that level are due to 
density variations, particulate scattering being assumed 
as negligible. 

Barrett and Ben-Dov (1967) obtained soundings ex- 
tending to less than 1 km, and assumed in some cases that 
attenuation was negligible, and in others that the ratio of 
the particle backs cat ter cross section to the attenuation 
cross section was a known constant, obtained from a model. 

o 

Furthermore, they assumed a molecular number density pro- 
file. 

With recent advances in laser technology and In the 

production of photomultipliers and optical filters, it has 

become possible to examine the spectral composition of the 
* 

light returned in a laser radar system, and thereby to ob- 
tain additional Information needed to dispense with some of 
the assumptions discussed above. 

The high monochromaticity of the laser source, used in 
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conjunctioi^with an optical filter with a sharp frequency 
response, enables the investigator to monitor regions of 
the spectrum extremely close to the frequency of the 
transmitted light, with minimal interference from that 
light. The high photon densities obtainable through use of 
lasers, together with the high quantum efficiencies affor- 
ded by modern photomultiplier tubes, allow the observation 
of secondary scattering phenomena that were hitherto un- 
observable outside of the laboratory. This secondary, or 
Raman scattering, comes about as a result of the exchange 
of energy '.between the incident photons and the matter com- . 
posing the atmosphere, and appears as light scattered at 
frequencies different from that of the incident light 
( Her zberg, 1939). The character of this exchange is 
uniquely dependent upon the molecular characteristics of 
the particular substance and, as such, can be used as a 
means of detecting the presence of that substance. 

Leonard (1967) was the first to demonstrate that I 
Raman scattering could be employed for atmospheric obser- 
vations when he monitored returns from atmospheric nitrogen 
and oxygen. These returns were the result of vibrational 

quantum transitions (vibrational Raman scattering), and 

-1 -1 

represented frequency shifts of 2328 cm and 1551 cm 
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due to and 0^ respectively. 

Nitrogen and oxygen also produce a rich spectrum of 
scattering frequencies due to rotational quantum transi- 
tions (rotational Raman scattering). These frequencies lie 

-1 -1 

in a region extending from about 10 cm to 200 cm to 
either side of the incident frequency. Salzman, et.al. 
(1971) are presently attempting to exploit the strong 
temperature sensitivity of the lines in the vicinity of 
200 cm*"* to obtain remote temperature prof ilesv 

In liquids and solids, where intermole cular forces 
are greater than in gases, molecular rotation is inhibited 
or absent entirely. In its place, however, there may exist 
molecular librations (in liquids) or lattice vibrations 
(in solids) which act to produce scattering in the same 
range of frequencies as molecular rotation. 

If it can be shown that the off-frequency component of 
light scattered by liquid and solid constituents of the 
atmosphere is negligibly small compared to that scattered 
by the gaseous components, and if a relationship can be 
found between the intensity of off-frequency scattering 
and molecular number density, then we should be able to 
remove at least part of the arbitrariness which now exists 
in the interpretation of density and turbidity profiles. 
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Schotland (1969) discussed this possibility in general 
terras* 

It will be shown here that the off-frequency component 
of light scattered by liquid and solid constituents of the 
atmosphere is, indeed, negligible compared to that from 
gaseous components. It will be shown that, through 
selective weighting of the received off-frequency compo- 
nents, a relationship may be established between these 
components and molecular number density. An investigation 
is conducted into the feasibility of developing a practical 
laser radar system which employs these principles to obtain 
density and turbidity profiles. 
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List of Symbols 

System constant for elastic channel (W~cra). 

System constant for Raman channel (W-cm). 

—3 

Molecular number density of air (cm ). 

Power received in elastic channel (W) » 

Power -received in Raman channel (W) • 

Atmospheric transmission - molecular component. 

Atmospheric transmission at frequency of elastic - 
scattering - molecular component. 

Atmospheric transmission at frequency of Raman 
scattering - molecular component. 

Atmospheric transmission - particulate component. 

Atmospheric transmission at frequency of elastic 
scattering - particulate component. 

Atmospheric transmission at frequency of Raman 
scattering - particulate component. 

Height (cm). 

/ _ ‘"‘1 V 

Molecular absorption coefficient (cm ). 

Attenuation coefficient - molecular component (cm *)• 

Attenuation coefficient * particulate component (cm ). 

particulate elastic backscatter coefficient of air 
(cm" 1 ). 
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*1 Particulate Raman backscatter coefficient of air 
r (cm" 1 ). 

Molecular elastic backscatter cross section of air 
e (cm^-ster“i) . 

2 -1 

Raman backscatter cross section of air (cm -s ter ). 
Atmospheric turbidity. 
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A» General, 

Consider a pulsed, monostatic laser radar system, with 
a receiver which simultaneously monitors the frequency of 
elastic scattering (i*e., the transmitted frequency), 
and a frequency V r * at which rotational Raman scattering 
occurs. As shown by Schotland (1969), the applicable 
radar transfer equations are 



Elterman (1.968) tabulated atmospheric attenuation 

profiles for the entire visible spectral region, using a 

large body of observed data, together with an attenuation 

model® Some figures derived from his results are presented 

in Table 1. The references to visibility are estimated 

following the discussion of meteorological visual range 

given by Middleton (1951), The table shows that, for 

visibilities greater than about seven miles, the variation 

of transmittance with frequency is generally below one 

•~1 

percent per 100 cm . For the very turbid atmosphere in 
Table l.c., the variation is still less than two percent 
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Table 1 . 

a, Variation with frequency of total one-way atmospheric 

the surface and selected altitudes (percent change per 1 00 cm ; 

for an atmosphere with attenuation coefficients equal XibrntfaTth* 
Elterman (1968). For this atmosphere, the horizontal visibility at th 

surface is approximately 1U miles* 

SKSfc it it it it it it It It it 


2 

(to) 

2 

5 

10 

20 

oo 


0,19 

0.33 

0.U7 

0.51 

0.61 


0.1U 

0.2U 

0.3U 

oil 

O.UO 


0.11 

0.09 

0.10 

0.09 

0.07 

0.07 

0,18 

0.15 

0.16 

0.11* 

0.10 

0.09 

0.30 

0.21 

0.21 

0.18 

0.13 

0.12 

0.26 

0.23 

0.31 

0.26 

0.16 

0.15 

0.23 

0.19 

o.UU 

0.3U 

0.18 

0.17 


0.09 

0.11 

0.13 

0.15 

0 . 1 $ 


b. Same as a., but with merman's attenuation coemci^ts^ltiplied by 

uniform factor of two (Surface visibility approximately 7 miles) . 


2 

.5 

10 

20 

OO 


0.38 
0.71* 
1 .09 
1.39 
1 .50 


0.27 

0.1*9 

0.73 

0.91 

1.15 


0.22 

0.38 

o.5o 

0.55 

0.75 


0.1? 0.21 0.22 0.23 0.27 

0.30 0.33 0.33 0.31* 0.39 

0 JJi 0.36 , 0.36 0 JiO Olio 

0.57 0.5V 0.55 0.56 0.62 

0.73. 0.70 0.69 0.75 0.68 


0.35 

0.1*5 

0.55 

0.68 

0.86 


_ s awe as a., but with Elterman 's attem&tion coefficients maltlplied by 

uniform" factor of five (Surface visibility approximately 2.8 miles). 

0.69 ' 0.88 

0.99 1 -18 

1 .19 1 -U3 

1 .55 1 -65 

1 .75 1 -70 


2 

5 

10 

20 

oo 


0.82 
2.12 
3.30 
I* .57 
5.51* 


0.73 
1 .3U 
2.08 
2.63 
2 .91* 


0.71 

0.U8 

0.53 

0.55 

0.57 

0.95 

0.82 

0.81 

0.79 

0.87 

1 UiO 

1.15 

1 .05 

1 .01 

1.16 

1.58 

1.51 

1 .1*5 

1 .1*8 

1.53 

2.19 

2.00 

1 .93 

1 .87 

. 1**1 
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for wavelengths longer than 0.60/f, and less than three per- 
cent per 100 cnf 1 for wavelengths longer than 0.45x/* It is 
evident that, given the range of altitudes and turbidities 

in which one wishes to operate, a range of wavelengths can 

* » 

be found for whi^hj 


(W ' *1 I 

I Jrae ■Jpe 


- 7 


me 


7„ e 7, 


pe *^mr ^pr 


% 


C 


= |7 7 ' 

1 •'m •'p 


(3) 


We shall assume for the present, and show in Chapter 
III.H.3, that 


n r « <r r n (4) 

« 

That is, the off-frequency component of scattering due to 
particulates is, at all altitudes, negligible in comparison 
to that due to molecules. 


B. Measurement of particulate backscatter. 

By combining (1) through (4) in such a manner as to 
eliminate we obtain 


% -**!• 

<T r n K e P r <T P 



the ratio l^/Kg is a system constant obtained through . 
laboratory measurements. P e and P r are the quantities 
measured by the laser radar system. The constant (T^il 
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is determined in Chapter IX I. G. and is an invariant quanti- 
ty for a given gas, as long as resonance scattering condi-v 
tions do not exist. 

Equation (5) is the expression which allows the de- 
termination of one of the quantities of interest here, 
fj / CT r n, which represents the ratio of light backs cattered 
by particulates to that backscattered by molecules, and is 
a measure of atmospheric turbidity. 

There is no standard definition of turbidity. For 
example, Elterman (1968) uses 

S 

r(z> - (6) 

1 d3L 

£p * Hz (7) 

„ 1 ■O'. . 

--r m *r <»> 

This definition includes absorbing properties of the con- 
« 

stituents as well as oblique scattering, while the quanti- 
ty/^ e / CT r n does not. A definition similar to that of 
Elterman (1968) is that of DeBary, et.al. (1965)* who used 


where 


and 
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1 + 


ln 


1 




The relationship between n and (J is simply 

3T w 



frs 

3 ®r 


<Tn + <* 
r m 


( 10 ) 


via Rayleigh’s phase function for small scatterers 
(Chandrasekhar, I960). In fact, there are many spectral 
regions accessible to laser radar for which 


8tr 

*- « 7 7 ** 

r 


(ID 

♦ 


making (10) even simpler. Unfortunately the relationship 
between and ^ is not as clear, depending as it does 
upon the particulate phase function and absorption coef- 
ficient, which are as highly variable as the size distrib- 
ution, shape and composition of the particles suspended in 
the atmosphere (Deirmehdjian, 1969). Thus, there can be no 
simple conversion between the ratio, <T r n, to be meas- 
ured here, and the various other definitions of turbidity. 

The primary requirement for a definition of turbidity 
to be useful from an operational standpoint is that the 
quantity can be measured reproducibly at different times 
and places. . The ratio 0^n satisfies this requirement. 
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C. Measurement of tropospheric transmission and particle 
attenuation* 

In equation (5) it was required only that we know the 
ratios K /K and (T/ffil* Suppose that we are able to obtain 
K^. and <F^ explicitly. Then, for each altitude z, using 
(2) through (4), we may solve the off -frequency transfer 
equation for the transmission: 

V p - 

Molecular number density, n, is the only unknown quantity 
on the right side of (12).. At tropospheric altitudes, it 
may be estimated to a fair degree of accuracy from statis- 
tical information or from recent synoptic data. Thus, 
using the off-frequency transfer equation together with an 
estimate of the molecular number density at a given height, 
we can compute the transmission from the surface to that 
height. It should be noted that any error in the estimate 
of n results in an error in the computed transmission that' 
Is only half as great. 

In a spectral region where molecular absorption is 
negligible, an estimate of n implies an estimate of 
and as may be ascertained from equations (8) and (10)* 


L K r C r n j 


( 12 ) 
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Without any further assumptions then, we compute *3^ from 
equation (12) and, consequently, ^p, using equation (7). 


p, Measurement of upper atmospheric number densities 
Several studies, such as those by Smith (1964), 

Nee (1964) and that in the Handbook of Physics and Space 
Environments (1965) have shown that there exists near a 
height of 7 km, an ”isopycnic level” in the atmosphere, 
where density has a minimum variation in time and horizon- 
tal space. Smith’s (1964) study shows, for instance, that 
at Buffalo, N. V., 99% of all density observations at 7 km 
lie within ±3.6 percent of the annual mean (e.g., a RMS 
deviation of *1.2 percent) for that station. This value is 
typical of those given by smith for other stations in 
temperate and polar regions. While no similar studies 
could be found for tropical regions, it is expected that 
densities there could be even more closely estimated from 
climatology, since sharp air mass contrasts do not exist 
there . 

Elterman (1968) shows that, at a wavelength of 0.7xf , 


Ct *3 ) 2 decreases by only about 117. between 7, km and the 
top of the atmosphere. About 60% of this quantity, or 6.6 
percent of the total decrease is due to aerosols. 



Writing (12) as 


16 


n 


W^mV 2 


( 13 ) 


we see that, if (^ C J ) 2 at 7 km is determined by the 

in p 

method of the previous section, using the annual mean value 
of n for the location in question, then we may determine n 
from J km to the top of the atmosphere to a fair degree of 
accuracy, even if the contribution to J ^ above that altitude 
is neglected altogether. We can improve the accuracy of our 
determination of n considerably by employing some model for 
^p, such as an extrapolation of ^(7 km) according to 
Elterman*s (1968) aerosol extinction profile. Elterman's 
aerosol transmission factors for 0.45^ and 0.70/{ are 
presented as fractions of the aerosol transmission factor 
at 7 km in Table 2. 



Table 2. Two-way aerosol transmission factors (0 - z) 
relative to the value at 7 km km)J • 

After Elterman (1968). 


z (km) 


(*)/ tan) 


%<«>/ J p (7ta») 


(0.45^0 


(0.70/d 


7 

8 

9 
#• . 

10 


1.000 

0.99a 

0.984 

0.978 


11 

12 


. 14 
. 15 
16 

17 

18 * 

19 

20 
21 
22 
23 


0.970 

0.965 

0.957 

0.951 

0.946 

0.940 

0.934 

0.929 

0.925 

0.921 

0.918 

0.916 

0.916 

0.914 
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Chapter III* 

CHARACTER OF THE ROTATIONAL RAMAN SPECTRUM OF AIR 
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a 



J 


» 


c 



d 
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-k 


ei 


' e X 


List of Symbols 

2 

Cross sectional area of sample cell (era ). 

Avogadro f s number (mole~*). 

Radial distance of atom from molecular axis of sym- 
metry (cm). 

Molecular rotational constant at vibrational quantum 
number, i (cm“^). 

Probability cf occurrence of transitions between 
quantum ctate j and j * . 

1 

Speed of light (cm-sec ). 

Molecular rotational constant at vibrational quantum 
number, i (ciq ). 

Diameter of tubing (cm). 

Partial pressure of H^O in the atmosphere (dyne-cnf^}. 
Gravity (cm-sec ). 

Statistical weight of rotational quantum state, j. 

Scale height (cm), 

Planck* s constant (erg-sec). 

Relative intensity of scattered light. 

Relative intensity of parallel polarized component 
of elastic scattering. 

Relative intensity of perpendicular polarized ’ 
component of elastic scattering. 
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I g Relative intensity of elastic scattering of both 
polarizations. 




i 

j 

r 

K 

k 

L 

M 


m 


m 


part 


N 

N o 

77 


n • 


Relative intensity of parallel polarized component 
of Raman scattering (0- and S-branches). 

Relative intensity of perpendicular polarized 
component of Raman scattering (0- and S-branches), 

Relative intensity of Raman scattering of both 
polarizations (0- and S-branches). 

Molecular*’ moment of inertia (gm-cm ), 

Vibrational quantum number of initial state. 

Rotational quantum number of initial state. 

Rotational quantum number of final state. 

Constant . 

4 

Boltzmann constant (erg-°K). 

Horizontal length of tubing (cm). 

Atomic mass (gra). 

Refractive index. 

Mass concentration., of particulate matter in the 
atmosphere (gm-cca )• 

—3 **. • 

particulate number density (cm ), 

—3 

Particulate number density at z *» 0 (cm ). 

Average number of particles viewed by spectrometer 
during- integrating time interval. 

Molecular number density (cm ) . 


n, Number density of molecules in rotational quantum 
J state, j (cixf^). 
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-3 

Number density of molecules with energy W (cm ). 

“2 

p Atmospheric pressure (dyne -cm )» 

Q Volume flow rate of air through spectrometer sample 

cell (cnr-sec*^) . 

Re Reynolds number* 

r Particle radius (cm). 

Maximum radius considered in particulate number 
density computation (era). 

r , Minimum radius considered in particulate number 
m density computation (cm). 

T Temperature (°K) . 

v Maximum speed of particle motion through sample 
cell (cm-sec** 1 ). 

v Mean speed of particle motion through ^spectrometer 

sample cell (cra-sec** 1 ). 

W(J) Amount of energy in a rotational or vibrational 
state of quantum number J (ergs). 

W^j) Amount of energy in a rotational state of quantum 
number j (ergs). 

W v (i) Amount of energy in a vibrational state of quantum 
number i (ergs). 

w Settling velocity of particles (cm-sec"* 1 ). 

z Height (cm) . 

*> 

U Mean molecular polarizability (cm J ). 

y Anisotropy of molecular polarizability (cm ). 

A Depolarization ratio. 

Frequency of Raman line relative to frequency of 
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elastic scattering (cnf*). 

At Spectrometer integrating time interval (sec). 

<Ty Spectral slit width of the spectrometer system (cm \ 

.Dynamic viscosity of air (cm -sec ). 

M Molecular mass of matter composing particulates 
(gra-mole - ^). 

V Frequency (cm ). 

<• 2 

Mass density of air (gm-cm ). 

—3 

^Part Mass density of particulate matter (gm-cm )* 

C T(7?) Standard deviation of the mean, 72, as viewed by the 

spectrometer. 

(T Raman backscatter cross section, including 0- and S- 
Razn branches and both parallel and perpendicular polari-’ 
nations (cm^-ster - * 3 *) . 

O' Rayleigh backscatter cross section, including Q- 

Ray branch Raman scatter and both parallel and per pen*? 
dicular polarizations (cm -ster~ i ). 
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A. General, 

in attempting to explain the scattering properties of a 
molecule, it is helpful to consider it as an asymmetric ag- 
gregate of electrons which act to deflect passing photons. 

The asymmetry of the electronic aggregate is a result of the 
asymmetry of bonding forces and the nuclear mass distribution, 
The molecular parameter usually used to represent this 
physical picture is the molecular polarizability, a tensor. 
For rotating molecules whose axes are randomly oriented 
with respect to the observer, it is convenient to refer to 
the polarizability In terms of its statistical properties: 

1. Mean polarizability, * - The portion of the 
polarizability that "remains constant as the molecule ro- 

tateSo 

2. Mean polarizability anisotropy, JC — The portion of 
the polarizability associated with molecular asymmetries. 

The deflection (or scattering) of photons due to the 
mean polarizability is thus independent of the shape of the 
molecule and is often known as Rayleigh scattering. We 
shall point out later, however, that the common operational 
definition of Rayleigh scattering also includes some por- 
tions of the scattering due to the polarizability aniso- 
tropy. Scattering due to the polarizability anisotropy 
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depends upon the asymmetries in the distribution of elec- 
trons, A molecule with perfect spherical symmetry has 
^ = 0. The intensity of rotational Raman scattering (that 

is, the scattering associated with molecular rotation) thus 
depends entirely upon the magnitude of if » 

In the sections which follow, we shall first discuss 
the characteristics of rotational Raman scattering, and 
then show more specifically how rotational Raman scattering 
and Rayleigh scattering are interrelated. 

The subject of rotational Raman spectra is one that is 
often fraught with many complexities. We are indeed for- 
tunate that the major components of the earth's atmosphere^ 
nitrogen and oxygen, exist overwhelmingly as homonuclear, 
diatomic, electronically unexcited molecules. Simple 
theoretical models of molecules such as these were among 
the early results of the development of quantum theory 
[[see, for instance, Placzek and Teller (1933)]. Numerous 
subsequent studies fe.g. Stoicheff (1959)3 showed Chat 
the simple models gave results that were compatible with 
actual observations. 

In recent years, with the increasing availability of 
sophisticated laser Raman spectrophotometers, second order 
Raman effects, not explained by the simple model have been 
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observed £e.g. A* 36 a nd Shimanouchi (1969)3 * These ef- 
fects constitute an extremely small fraction of the total 
Raman energy and have been observable only in conjunction 
with the rotational Raman lines that are very close to the 
exciting line. As we shall see later, these first few 
Raman lines are of only peripheral importance to our study, 
and second order effects superimposed upon these lines are 
of no concern to us. 

The theoretical development that follows will utilize 
the simple ••dumbbell 11 model for diatomic molecules. An 
extensive discussion of rotational Raman scattering in 
general, and this model in particular, is given by 
Herzberg (1939).* 

B» Population of molecular energy levels. 

At any given temperature, molecules of gases take on 
distributions of discrete rotational, vibrational and 
electronic states. Since we are concerned here only with 
neutral gases and" relatively low temperatures (200 K-320 £) , 
we shall consider only the ■••ground 1 * electronic state* 

The number of molecules in any rotational or vibration- 
al state depends upon the following factors: 

1. Boltzmann factor. If W(J) is the amount of energy 
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in a rotational or vibrational state of quantum number, J, 
then according to the Boltzmann energy distribution law 
(Herzberg, 1939), the number of molecules, dn^, that have 
a rotational or vibrational energy between W and W+dW is 
proportional to e "’W(J)/kT < j^ We see that, as temperature 
is increased, more molecules are found in higher energy 
states, relative to the number found in lower states. 

2. Degeneracy. Molecular motions are ” space quan- 
tized." That is, only certain modes of motion are permit- 
ted, with each mode having an equal likelihood of occur- 
rence. It can be shown (Herzberg, 1939) that 2J+1 times 

f-V. 

as many modes are permitted for a molecule in the J 1 state 
of vibration or rotation as in the zeroth state. 

3. Statistical weight, g^. Due to the influence of 
nuclear spin within molecules upon the selection rules for 
the occurrence of the various rotational quantum states, a 
regular alternation in the relative number of molecules In 
each state takes place (Herzberg, 1939). In homonuclear, . 
diatomic molecules such as oxygen, in which the nuclei have 
no spin, the ratio of the number of molecules in even- 
numbered states to that in odd-numbered states. is 0:1. 

That is, even-numbered states are completely unpopulated. 

In the case of nitrogen, where nuclei have a spin, the 
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ratio is 2:1. The even-numbered rotational states of ni- 
trogen have twice the population of odd-numbered states. 

Combining the above factors, we may state that the 
ratio of the number of molecules in the i t ^ 1 vibrational 
state and rotational state, n. to the total number’ 
of molecules, n, is: ' 



The sum in the denominator of (14) represents summation 
over all possible vibrational and rotational states. 
Herzberg (1939) showed that we may write, to a good 
approximation, 

W R = [ B i + X > " D i j 2 <j + U 2 (15) 

where and are molecular rotational constants which 
depend somewhat on the ' vibrational state, 1, of the mole- 
cule. Since the discussion here will concern itself only 
with the pure rotational Raman spectrum (i * 0), we shall 
concern ourselves only with 

W R = [ B o + 1) “ 0 o j 2 (J + l) 2 



(16) 
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The magnitude of W R is such that, at atmospheric tempera- 
tures, rotational states in the range, 0<j<30, are sig- 
nificantly populated, Values of B and D for N« and CL 

oo22 

have been computed from laboratory observations by 
Stoicheff (1959) and are known to sufficient accuracy for 
the present application. 

Concerning the character of W v , we shall state only 
that it is much larger than W R , so that at atmospheric 
temperatures, the ratio of the number of molecules in the 
first to that in the zeroth vibrational state is about 
l,40xl0~ 5 for and 5 „ 74 x 10"^ for 0^ 9 as observed by 
Herzberg (1939), That is, virtually all atmospheric l* 2 
and <>2 molecules are in the zeroth vibrational state. 
Thus, (14) can be simplified to read 



i=0 (17) 


in the atmospheric range of temperatures, 

C, Rotational Raman transitions. 

The rotational Raman effect is the result of transi- 
tions of the type. 



j - 


0 ,,- £2 
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where j is the quantum number of the initial state and j 1 
that of the final state. 

When the collision between a photon and a gas molecule 
is elastic, then a transition of the type, j f - j = 0, 
Occurs. The Raman- scattered radiation is thus at the same 
frequency 0.0 trident light. The lines of scattered 
energy (one for each j-valuo) r e . su i ting f rom this type of 
transition are superimposed at the frequency o* c he Inci- 
dent light -and are known as the ,T Q~branch M . They are some- 
times described as that portion of the Rayleigh scattering 
due to the anisotropy of the polarizability tensor. 

When a collision is accompanied by a net transfer of 
energy to the molecule, then a transition of the type, 

j ! - j = +2 

occurs, and the Raman-scattered radiation is lower in 
frequency than the incident light. The lines in which 
energy is scattered at frequencies lower than that of the 
incident light are known as "Stokes'* or "S-branch" lines 0 

For the case where there is a net loss of energy by the 
molecule, we have, 

j 1 - j = -2 

and the Raman scattering occurs at frequencies higher than 
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chat of the incident light. These lines are known as 
"anti-Stokes" or "O-branch" lines* 

The frequency separation between successive lines can 
be closely determined by considering the molecule as a 
"rigid rotator"; that is, by assuming that the centers of 
mass composing the molecule remain fixed relative to each 
other as the molecule rotates. When this assumption is 
made, the lines are found to be separated from the frequen- 
cy of the incident light by an amount 

AV = ±4B 0 (j + 3/2) (18) 

as shown by Herzberg (1939). For large values of j (large 
rotation rates), the "rigid rotator" assumption becomes 
less accurate, and a correction term must be added to (18) 
to account for the effect of centrifugal distortion. The 
expression then usually takes the form, 

AV = ± J^(4B 0 - 6D q )( j + 3/2) - 8D 0 (j + 3/2) 3 J (19) 

Neglect of the correction term in (19) results in less than 
a one percent error in Av at j~30, for both nitrogen and 
oxygen. 

The intensity of the j th rotational Raman line is shown 
by Herzberg (1939) to be proportional to the Boltzmann fac- 
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Cor and the statistical weight of the initial level 


g e " W R (j)/kr 

and the degeneracy of the lower of the two levels, 


(2j + 1) 
Stokes 


or 


(2 j* + 1) 
anti-Stokes 


as well as the transition probability, b^, , and the square 
of the anisotropy of the molecular polarizability, * 


That is. 




(2j + 1) gj bj, e-%<j> /kT y2 

_ TWi + 1) Si bt e-^ 57 ^] 11 

0,Q,S i’ol J J 


Stokes 


and 




(2 j ' + 1) gj bj, e"% (j)/la 2 

X + O &] b t e“V j)/kT | n 


anti-Stokes 


( 20 ) 


( 21 ) 


Transition probabilities, b~! , , have been derived by 
Placzek and Teller (1933), The expressions corresponding 
to pure rotational Raman scattering are: 
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J _ 3(„ i + , ... 

b j> 2(2 j + l)(2j +3) 

(S- and 0- branches) 


( 22 ) 


and 



1U +A) c 

(2j - 1) ( 2j +3} 


(Q- branch) 


(23) 


Values of X for many substances have been determined 
in the laboratory, primarily through measurements of de- 
polarization ratios. Such measurements were made by 
Bridge and Buckingham (1966) for and 0 2 * Their results 
were checked by this author through comparison of their 
figures with measurements of the rotational Raman spectrum 
of air made by Miller, et.al. (1968). If the ratio of 
molecular number densities of N 2 and 0^ is known, then the 

ratio can be extracted from the observed spec-* 

n 2 °2 

trum of air through comparison of the intensity of *? 2 
lines relative to the 0^ lines. This was done, and the 
ratio y ^ 2 , so determined, differed from that ob- 

tained by Bridge^and Buckingham by only a fraction of a 
percent. 

If we neglect, for the time being, the dependence 
inherent in all the intensities (Goody, 1964), then sever- 
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al statements may be made about the distribution of energy 
scattered in the pure rotational Raman spectrum: 

1. At atmospheric temperatures, where vibrational 
transitions are insignificant, the total energy scattered 
in the 0-, Q~ and S-branches ramains constant with changes 

QO 

in temperature. This may be verified by computing jL Ij 
for the 0*”, Q- and S-branches, using Equations (20) through 

(23). 

2. In a like manner, we find that the fraction of the 
energy scattered in the Q-branch remains constant at one 
fourth of the total. Thus, the total energy scattered in 
the 0- and S-branches must be independent of temperature. 
When temperature is .increased, scattered energy shifts from 
lower to higher j-values, and from the S-braneh to the 0- 
branch. 

3. The S-branch is always more intense than the 0- 
branch, since scattering intensity is proportional to the 
number of molecules in the initial state, and there are al- 
ways more molecules in the lower- lying qauanturo states than 
in the higher ones (see Section III.B.l.). From computa- 
tions to be made later in this chapter, we shall see that 
for air, the average S/0 intensity ranges from about 1*50:1 
at 200°K to 1.39:1 at 300°K. 



33 


D„ Raman scattering in relation to Rayleigh scattering. 

Having determined the intensities of the various rota- 
tional Raman lines relative to each other, we wish to ex- 
amine the rotational Raman scattering in relation to 
Rayleigh scattering. This relationship is best approached 
from a classical point of view, although it must always be 
recognized that classical mechanics does not fully explain 
the Raman scattering phenomenon. The classical derivation, 
for instance, predicts a continuum of scattered energy, 
rather than a series of discrete lines. Furthermore, it 
predicts that the 0- and S-branches will contain equal 
amounts of energy. Never the less, this technique does pre- 
dict the total energy in the 0- and S-branches relative to 
that found in the Q- branch, and that found in the portion 
of the Rayleigh line due to the isotropic part of the 
polarizability tensor, 

Stuart (1934) performed the computations in detail* 

For 90-degree scattering, in the plane perpendicular to 
the plane of oscillation of the electric vector of the in- 
cident field, the relative intensities, per unit solid 
angle , are z 
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+ 7y 2 /180 


(24a) 

(24) 

(24b) 



V 6 


X 2 / 15 

y 2 /2 o 



(25a) 

(25) 

(25b) 


It can be. shown that the same results hold for 180-degree 
scattering in that plane. For a direction making an angle, 
0, with the parallel direction, the intensity is simply 
given by. 

I Q - l tl cos 2 e + Ij_ sin 2 e 

as is shown by' Chandrasekhar (I960), We thus have a com- 
plete picture of the relative intensity and angular depend-* 
ence of rotational Raman and Rayleigh scattering. This is 
shown pictorial ly in Figure 1, 

Since we know the distribution of energy within the 0- 
and S-branches from earlier calculations, we are now in a 
position to compute the intensity of any individual rota- 
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2 



Figure 1, Angular dependence of rotational Raman and 
Rayleigh scattering® 
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tional Raman line relative to that of the Rayleigh line, for 
any scattering angle. 

Further, we may compute the total scatter cross section 
for any line (that is, the cross section integrated over 
all solid angles). 

E. Rotational constants and computations for nitrogen. 

Stoicheff (1959) computed the rotational constants for 
from observations of rotational Raman line spacing. His 
values are: 

B q = 1.9897 ± 0.0003 cm” 1 
D q = (6.1 ± 0.5) x 10” 6 cm” 1 

He points out that these values are in excellent agreement 

with values computed from the electronic spectrum by 

Wilkinson (1957), namely 

B = 1.9898 ± 0.0005 cm” 1 
o 

D « 6.4 x 10 ~k cm~* 

o 

Stoicheff 's' values are used in the present calculations. 
Relative line locations and intensities within the 0~, Q- 
and S-branches of the nitrogen spectrum computed using 
Equations (19). through (23) are presented in Figure 2. 

The relative intensities calculated here are compared 
to those observed by Miller, et.al, (1968) in Figure 3‘. 
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Figure 2. Computed rotational Raman line locations and intensities for nitrogen.- 
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Miller, et.al., in their observations, made no attempt to 
determine the gas temperature. Thus, in order to compare 
the present computations with their observations, it was 
first necessary to estimate the temperature at which the 
measurements were made. To accomplish this, line inten- 
sities were computed at a number of temperatures, and the 
results for each line plotted against observed values ex- 
tracted from chart records given by Miller, et.al* For 
incorrect temperatures, the result was a graph similar to 
Figure 3, but with the points representing high- and low- 
numbered lines scattered to opposite sides of the best fit , 
straight line. The temperature value inserted into the 
computations was adjusted until the selective scatter of 
points was minimized. The final result, shown in Figure 3, 
is a series of points falling close to a straight line 
distribution. A straight line relationship between the two 
sets of intensities shows that the observed and computed 
intensities differ from each other by a constant factor 
(that is, the shape of the two distributions is the same). 

F. Rotational constants and computations for oxygen.. 

Table 3 gives rotational Raman line spacings observed 
in 0 ? by Stoicheff (1959), Through use of the graphical 
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Table 3. Rotational line spacings in 0 2 determined by 
Stoicheff (1959). 


j 

AV (cm"* 1 ) 

1 

14,31 

3 

25,85 

5 

37.37 

7 

48.84 

9 

60.33 

11 

71.80 

13 

83.26 

15 

94,72 

17 

106.16 

19 

117,57 
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technique described by him, the applicable rotational 

constants were determined by this author. They are: 

B = 1.439 cm -1 
o 

D = 4.05 x 10 -6 cm" 1 
o 

For purposes of * comparison, B q was also calculated using 

independently computed constants: the radial distance of 

the oxygen atoms from the molecular axis of symmetry, 

-9 

a * 6,037 x 10 cm 

which is given by Wacker, et.al. (1964), and the mass of an 
oxygen atom, 

M = 2.656 x 10" 23 gra 

This allows us to compute the molecular moment of inertia, 

= 2Ma^ ~ 1,936 x 10 39 gm era 2 (26) 


From classical mechanics, we have the expression for the 
rotational constant, 



(27) 


Inserting our computed value of ^ 


we obtain 


B q — 1,445 cm 

which differs from the value computed from Stoicheff f s 
(1959) data by about 0,4 percent. 
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the values of B Q and D q computed from Stoicheff's 
data are used here. Relative line intensities computed 
using (19) through (23) and line locations within the 0-, 

Q- and S-branches of the oxygen spectrum are presented in 
Figure 4. The relative intensities calculated here are 
compared to those observed by Miller, et.al. (1968) in 
Figure 5. A detailed explanation as to how this comparison 
was carried out is given in the previous section, in the 
discussion of Figure 3. 

G. The rotational Raman backseat ter cross section of 
air (N„ and 

Through use .of the expression derived by Chandrasekhar 
(1960) for the Rayleigh backscatter cross section, together 
with the results of Stuart (1934) described in Chapter 
III.C, for the ratio of Raman to Rayleigh intensities, we 
determine that the rotational Raman backscatter cross sec- 
tion of air is given by 


®Vam^ lr) “ (t) [^ y2(Air) . 


(28) 


< 3 > + S 

The effective polarizability anisotropy of air, squared, 
3f 2 (Air), is merely an average of squares of the aniso- 


tropies for and 0^, weighted by number density: 
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FREQUENCY RELATIVE TO TRANSMITTED FREQUENCY (cm-') 

Figure 4. Computed rotational Raman line locations and intensities tor oxygen. ■ 
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JT 2 (Air) 


ir 2 (» 2 > 


x—s 

CM 

O 

CM 

n(N 2 ) 

y z (N 2 ) 

n(0 z ) 


n(N 7 ) 

1 

+ — 


n(° 2 )^ 


(29) 


The number density ratio in the atmosphere is cons id- 
ered in the U« S. Standard Atmosphere, 1962 to be constant 
with height between the surface and 90 km* The number 
density ratio is given as 

n(tt 9 ) . 78.084 

•« =- « - 3.728 (0 4 km) 

n(0 2 ) 20.946 


Values for the polarizability anisotropies in and 

0 2 were obtained by Bridge and Buckingham (1966) from 

measurements of depolarization ratios. Their values, which 

o 

were obtained from measurements made at 6328 A, are 
tf(N 2 ) - 0.69 A 3 
<0 2 ) = 1.09 A 3 

These figures give a ratio. 


V 2 (o 2 ) 


~ 0.401 


The Landolt-Bornstei n Tables (1951), citing observations 

mu aoaac i j i uwm*iniw | » , » l 1 M " 


made in 1928, give values of 
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V(N 2 ) - 0.93 

y<0 2 ) = 1.14 A 3 


which result in the ratio, 


% 2 {n 2 ) 


0.665 


Ory and Yura (1965) computed # for N 2 from theoretical 

°3 

considerations and obtained a value of 1.02 A . 

There are obviously, considerable differences among the 
various determinations of polarizability anisotropies. 
Although the determination of quantitative turbidity sound- 
ings depends upon knowledge of the absolute values of these 

anisotropies, the design considerations which follow depend 

2 2 

only upon the ratio 2 (N 2 )/}f*(0 2 ). 

Miller, et.al. (1968) observed the composite rotational 

Raman spectrum of air directly. Through the use of the 

relative strengths of N 2 and 0^ lines in their results, to* 
• • 

gether with the precisely known ratio, n(N 2 )/n(0 2 ), the 
ratio y V ^(0 2 > was computed here for comparison with 

the above figures. 

The observations by Miller, et.al. (1968) were found to 
agree closely with those of Bridge and Buckingham (1966) 
(See Appendix 1.), and thus their figures are used in the 
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present study. We therefore find that 

y 2 (Air) = 1.317 y 2 (N 2 ) 

and, for Bridge and Buckingham’s value of Y (N^) 3 

Y 2 (Air) = 0.627 A 6 


This gives the rotational Raman backscatter cross section 
of air as 

0*1 (Air) — 4,91 x 1 cm 2 -ster ^ (A- 6943$) 

Ram 

8f*X 

o+g 

At 30G°K, the distribution of scattering between the S~ and 
0-branches is in the ratio 1.39:1, or 

2 —1 

Cr n (Air) - 2.85 x 10 cm -ster 

Ram 
K+X 

and . S 


0^ atD (Air) ~ 2.06 x 10 30 cm 2 - star’ 1 

jf+x 

© 

Figure 6 shows the composite rotational Raman spectrum of 
air. 


H. The Rayleigh backscatter cross section of air (N „ 

The Rayleigh backscatter cross section of air has been 
variously defined to include all of the rotational Raman 
lines, to include some of the rotational Raman lines, to 
include only scattered radiation polarized parallel to the 
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Incident light, or to include both parallel and perpen- 
dicularly polarized radiation. 

Each of the above definitions has merit and has enjoyed 
some popularity in the literature, and it is of course 
necessary to keep these various definitions in mind when 
comparing the results of different investigators. The 
definition of the Rayleigh backscatter cross section to be 
used here includes only light returned at the wavelength 
of the incident light, without regard to polarization. As 
such, through use of the expression derived by Chandrasek- 
har (1960) together with the results of Stuart (1934), we 

’ * 

find the Rayleigh backscatter cross section to be 


*1 Ray< Air) 

■/Hi. 



cr(Air) 


+ 2 ( Air ) 


(30) 


In a manner analogous to that of the previous section, we 
obtain 


<* 2 (Air) - C< (K z ) 


c< 2 <0 2 ) n(K 2 ) 


e< 2 (N 2 ) 


n(0 ) 


n(N 2 ) 

n(0 2 ) 


(31) 


The polarizability of air may be determined from 
measurements of the refractive index and total molecular 



50 


number density by means of the Lorentz - Lorenz formula 
(Born and Wolf, 1959), 

o< = 

Elterman (196B), using an expression determined by 

Edlen (1953), gives the refractive index of air at 6943 8* 

for standard temperature and pressure (1013 mb and 15°C) as 

m *= 1.00027584 

The U. S. Standard Atmosphere, 1962 gives the number densi- 
ty at that temperature and pressure as 

19 -3 

n « 2.5471 x 10 * cm J 

Thus , we obtain . 

*7 ofi 

, <X' z (Air) = 2.97 A 

from considerations of refractive index. 

The Landolt Boms te in Tables (1951) give 

C<(N 2 ) = 1.76 A 3 

o<(0 2 ) = 1.60 A 3 

from spectroscopic measurements. These figures result in 

7 °6 

CV (Air) - 2.99 A° 

Stans bury j Crawford and Welsh (1953) computed OCO^) 

and 0< (0^) using refractive indices given in the Inter - 

o ' 

national Critical Tables for 4358 A. They give 


m 2 - 1 
m 2 + 2 


(32) 
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?3 


C<(N 2 ) - 1.79 A 


c<( o 2 ) 


O'* 

1.63 A J 


or 


<* 2 (Air) = 3.10 A 6 

Ory and Yura (1965) derived o((N 2 ) from molecular orbital 
considerations. Their theoretical approach yielded 

<*(N 2 ) = 1.41 a 3 

Kolker and ‘Karp lus (1963), who also employed a theoretical 
approach, found that 


<*(n 2 > 


2.40 A 3 


In this study, we adopt a value for (Air) intermediate 

between that computed from refractive index considerations 

06 . 

(2.97 A ) and that computed from spectroscopic data 

06 

(2.99 A ), or 


c( 2 ( Air) 


2.98 




It is felt that this value is representative of our present 
knowledge of c<^(Air) . 

2 

Using this value of C< (Air), together with the value 

r\ 

of X (Air) found in the previous section, we obtain 

„0 Q 2 

(Air) =: 2.01 x 10 cm ~ster 

#I+£ 

o 

for a wavelength of 6943 A. This value is 40.9 times as 
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large as the Raman backscattef cross section of air. 

I. Of f-frequettcy components of scattering from minor 
constituents of the atmosphere. 

1. Gases. 

a. Permanent constituents. 

Table 4 shows the relative abundance of gaseous atmos- 
pheric constituents with molecular number densities great- 
er than 10 of the total. The monatomic gases. A, Ne, He 
and Kr, by virtue of their inherent molecular symmetry, have 
no rotational Raman spectrum (Herzbcrg, 1945). Methane, 

CH^, while not monatomic, also has no rotational Raman 
spectrum, since it is a "spherical top" molecule, and 
therefore has no polarizability anisotropy (Herzberg, 1945). 

Thus, the only permanent gaseous atmospheric constitu- 
ent other than and 0^ which could produce rotational 
Raman lines is C0 2 , a linear triatoraic molecule which is 

highly anisotropic. Bridge and Buckingham (1966) give the 

°3 

polarizability anisotropy of C0 2 as 2*10 A . 

We may write the ratio of the intensity of Raman scat- 
tering from CO ^ to that from and 0^ as 

I(C0,) n(C0 9 ) iT 2 (C0„) 

■ = — — — s (33) 

I( V°2> n(N 2 )Jr‘ i (N 2 ) + n(0 2 )^ Z (0 2 ) 



Table 4. Relative abundance of permanent gaseous 
constituents of the atmosphere with molecular 
number densities greater than 10 ^ of the total, 
(After Goody, 1964) 


Molecule 

Fraction of Total 
Number Density 

N 2 

7.8084 x 10" 1 

° 2 

' 2.0946 x 10 -1 

• A . 

9.34 x 10" 3 


-4 

C0 2 

3.3 x 10 

Ne 

1.81.8 x 10" 5 

He 

5.24 x. 10*' 6 

CH, 

4 

1.6 x 10“ 6 

Kr 

1.14 x 10“ 6 
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When we insert the appropriate anisotropies and number, 
densities into (33), we arrive at an intensity ratio of 
0*23 %. Thus, for the purposes of this study, we may neg- 
lect the effect of rotational Raman scattering from CO^, 
as well as that from the other permanent constituents of 
air other than N 2 and C> 2 » 

b. Water vapor* 

(1) Rotational Raman backscatter cross section. 

H 2 0 is a planar asymmetric top molecule and, 
as such, has a very complex rotational Raman spectrum 
(Kerz.berg, 1945)* We shall not concern ourselves here with 
individual spectral lines, but only with the ratio of the 
total water vapor Raman cross section to that of air. How- 
ever, it should be noted here that Eisenberg and Kauzmann 

(1969) state that, at room temperature, the rotation band 

~1 

is most intense near AV = 200 m * 

In Chapter III.F., we determined the polarizability 

anisotropy, , of air, and noted that rotational Raman 

scatter cross sections are proportional to if independent 

of the nature of the spectral distribution of the lines or 

y 2 

molecular geometry. We thus limit our task to finding 6 
for in order that we may compare it to if 2 for air* 
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The optical properties of water vapor are not well 
known* As a simplification, Dorsey (1940) assumed that 
two of the axes of the polarizability ellipsoid are equal. 
By employing this approximation, he was able to treat the 
H 2 0 molecule as though it had an axis of symmetry (that is, 
that the polarizability ellipsoid was an ellipsoid of 
revolution) . For such a molecule, the polarizability an- 
isotropy is given by Stuart (1934) as 



(34) 


°3 

Dorsey (1940) gives o( as 1.59 A and Fabelinskii (1968) 
gives A as 0*. 020. Thus, we find that . 


Jf 2 (H 2 0) - 


°6 

0.39 A 


2 v 6 

In Chapter III.G«, it was shown that if (Air) is 0.629 A . 
Then the ratio of the squares of the polarizability aniso- 
tropies for H 2 0 and for air is 0.62. 


(2) .Number density in the atmosphere. 

In order ‘to place limits on the water vapor number 
densities which might be encountered in the atmosphere, we 
recognize that atmospheric water vapor content decreases 
much more rapidly with height than does the number density 
of air. As such, we expect to find the worst case at or 
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near the surface. Consider a very moist atmosphere with 
p = 1013 mb, T *■ 30°C and RH « 100%. Then, using the 
table given by Hess (1959), we obtain 

n(R,0) e 4.243 x 10 4 

£ — = _ « = 3.99 x 10 * (35) 

n(Air) p 1.013 x 10 

or 3.99 percent. 

(3) Raman scattering from water vapor relative to that 
from air. 

The ratio of scattered light intensity from H 2 0, rela- 
tive to that from air is given by 


ICH 2 0) nO^O^CHgO) 

I (Air) n(Air) i "(Air) 


(36) 


Using the results of Sections (1) and (2) above, we find 
that 

kh °) ; m2 

* — - 2.47 x 10 

I(Air) 

or 2.47 percent, under the pessimistic circumstances 
cited above. 
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2 . Liquids ♦ 
a* General. 

When a gaseous substance condenses into liquid form, 
intermolecular forces become much greater and, as a result, 
molecular rotation is hindered or absent entirely. 

Herzberg (1945) discusses the nattire of Raman scattering in 
liquids. He shows that, instead of a series of lines, the 
off-frequency component consists of a continuum of energy 
extending to either side of the undisplaced line, and de- 
creasing monotonically with AV* Superimposed on this 
continuum may be broad lines resulting from transitions 
between hindered rotational or translational states and 
intermolecular bending and stretching modes. 

The only liquid present to any appreciable extent in 
the atmosphere is water, and the remainder of this dis- 
cussion will be confined to this substance. 

b. Raman spectrum of liquid water. 

The Raman spectrum of liquid water has been studied by 
several investigators [e.g,, Gross*, et.al, (1937), 
Narayanaswamy (1948), Walrafen (1964, 1967) and Senior and 
Thompson (1965)J , It is generally agreed that, % a ad- 
dition to the continuum mentioned in Part a. above, broad 
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*■1 ■» 

lines exist near = 60 cm* and AV — 165 cm > as well 
as at other frequencies outside out region of interest* 
Unfortunately, none of the above studies deals with the 
intensities of these lines relative to the elastic scat- 
tering* The study by Walrafen (1964) shows that both lines 
decrease in intensity by. about a factor of three as tempera- 
ture is increased from 0°C to 40°C. 

To place an upper limit on the Raman scatter cross sec- 
tion per molecule of liquid water, we recall that such 
scattering is merely a hindered form of the scattering that 


occurs in water vapor and, .as such must be lower in inten- 


sity* Thus, we may write, as an upper limit, that 


cr - gi 

u Ram Ram 

H^O Liquid . H^O Vapor 


(37) 


Since we found previously that 


^Ram 
H 2 0 Vapor 


0,62 < r x 


we may write 


^Rara 

H2O Liquid 


o-62 cr; 


(38) 


(39) 
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c. Molecular number density of liquid water in the 
atmosphere. 

Typical clouds have water contents of the order of 

3 x 10 7 gm-cnr 3 , or about 1 x 10^ molecule s -cm 3 '. Air 

' * 19 

density ranges from about 2,5 x 10 molecules at the sur- 

1 R —3 

face to 8*6 x 10 molecule s -cm - at 10 km. Thus, a con- 

servative upper limit for the ratio of liquid water mole- 
cules to air molecules ranges from 4‘ x 10* at the surface 
to 1*2 x 10“ 3 at 10 km. Above 10 km, water clouds should 
be of no significance. 

Thus, for a relatively dense cloud at 10 km, under 
these pessimistic assumptions, the scattered intensity 
ratio would be 

Liquid t °H20 Liquid ^ 

( T v n(Air) 

(0.62) (1.2 x 10" 3 ) = 0.074 % 


When a substance assumes a solid state, molecules are 
bound tightly together into a relatively fixed lattice. 
Molecular rotation is not possible, but low frequency lat- 


^Cloud 

*Air 

3. Solids. 
a. General . 
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tice vibrations occur which cause off-frequency scattering 
components in the same region as rotational Raraan scatter-: 
ing* For a given substance, the character (frequency and 
intensity) of the lattice vibrations and, therefore, of the 
scattering, vary with lattice defects and impurities. As 
in the case of liquids, these lines are generally quite 
broad. Mitra (1969) presents an extensive discussion of 
Raman scattering in solids. 

b. Particulates in the atmosphere. 

Consider a highly turbid urban atmosphere in which the 
particle number density distribution varies with height 
according to 

z 

ft (particles-cm * ) « e (41> 


where the particle number density near the surface, is 
about 5 x 10 particles-cm and the scale height, H, is 
1.67 km. This is similar to several distributions shown by 
Junge (1963). We also assume that the size distribution is 
of the form. 


= Kr ~ 3 

d log r 


(42) 


which is given by Junge (1963). In alternative form, this 
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may be written as 


dN 

dr 


-4 
3 Nr ■ 

r min ■ ™ r max 


( 43 ) 


If we assume that the particles are spherical, then the 
number of molecules per increment of radius is 


di w 


part a dN 


(44) 


dr 


M 


dr 


Combining (43) and (44) and integrating, we obtain the num- 
ber of molecules found in particulates in the finite inter* 


val of radii, r min < r i 


max* 


^art 


4<rr W aN ln ;w 


*^ r min ~ r max * 


(45) 


min 


Let r m . = 10“ 6 cm and r “ 10 ^ cm, M - 100 gm-mole 
min max 

-3 


and ^ part =* 1’gm-cm* . Then, 


♦^art 


= 3.5 x 10 N 


5 -3 

At 2 » 0, where N = 5 x 10 particles^cm , 

11 «3 

np art « 1.73 x 10 molecules -cm 


The mass loading associated with- a particulate distribu- 
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tion of this sort is 

•^art 

" ,part ~ a 


28.8 vWg-m -3 (46) 


The distribution of mass loading with particle radius in- 
terval is shown in Figure 7. These quantities are of the 
same magnitudes as those found in urban areas. 

Distribution of particulate concentration with radius 
interval is shown in Figure 8. -Particulate concentrations 

..4 

computed in this manner for the region r>0.3 xA will be 
shown in Chapter III.H.3.c.(2)(a) to agree with observed 
concentrations . 

The molecular number density of air at z = 0 is typi- 
cally 2.5 x 10*^ molecule s-cm*^, as was mentioned earlier * 
Thus, for off-frequency scattering from particulates to be 
even one percent of that from air molecules, the scatter 
cross section per molecule would have to be 1.5 x 10** times 
the rotational Raman cross section of air. Were this the 
case, Raman scattering from large crystals would be an 
easily observable phenomenon, which it is not. Mitra 
(1969), in fact 9 states that Raman scattering from crystals 
is a second order process. 

By virtue of the molecular number densities involved, 
and the fact that Raman scattering from crystals is a 



PERCENTAGE OF PARTICULATE MASS ABOVE INDICATED R£ 



Figure 7. Distribution of particulate mass loading with particle radius interval, 

as eiven bv model discussed in tevt* . 


NUMBER OF PARTICLES WITH RADII ABOVE THAT INDICATED 
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second order process, we expect to be justified in assuming 
that the of ^frequency component of scattering from atmos- 
pheric particulates is negligible in comparison to that 
from air. In the next section, an experiment will be de- 
scribed which shows that this is indeed the case. 

c. Verification in the laboratory that the contribution of 
off-frequency particulate scatter to the Raman spectrum of 
air is insignificant compared to that of molecular scatter. 
(1) Design of the experiment. 

We wish to duplicate in the laboratory, using a Cary* 
Model 82 Raman Spectrophotometer, the ratio of off-frequen- 
cy particulate to molecular scattering observed in the at- 
mosphere, by a laser radar system. According to Haber (1972), 
the sample volume of the Model 82 may be considered as a 
cube, 50// on a side (l.e., a volume of 1.25 x 10~ 7 cm 3 ). 

If we assume that a typical urban atmosphere may contain as 
little as 10 particles per cubic centimeter, then, at any 
given time, there will be an average of only 0.00125 par- 
ticles in the sample vbltnne. Thus, we cannot rely on spa- 
tial averaging of particle number density to give us a 


*Cary Instruments, 2724 3. Peck Rd., Monrovia, CA 91016. 
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statistically significant result. We must instead rely 
upon time averaging as a means of obtaining statistical 

stability. 

Suppose that particles are forced through the bottom 
boundary of the sample volume with velocity v. Then the 
average rate at which particles pass through that boundary 

is 

*2L = KAv (47) 

dt 

where N is the average particulate number density, and A. 
is the area of the bottom of the sample volume. 

The average number of particles seen by the spectro- 
photometer in the integrating time interval is 

72 - - nav St (48) 

dt 

If we assume Poisson statistics, then the standard devia-^ 
tion of this number is 

&Oi) =7 i* = [NAvSt]^ (49) 

and 

jrsni = r N Av5ti-* (so) 

7t L J 

The integrating time, is. related to the spectral slit 
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width, S V , by 


St - 


dv /dt 


(51) 


where dv /dt is the spectral scan rate of the instrument 
(constant) . Thus , 

% 


fin) = f d^ /dt_T 
92 |_NAv Sv J 

Let £*> - 2 cm" 1 , d* /dt - 0.4 cm^-sec* 1 , and 

A - 2.5 x 10~ 5 cm 2 . Then 


*S 2L l 

72 


8 x 10 


41% 


L Nv 


or 


8 x 10 


N [<r(7l)/^] 

If N = lu* particles -cm , and 
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o.i 


Then 


«t 


v 


(52) 


(53) 


(54) 


800 cm-^sec 

If the sample cell through which this flow takes place has 
sectional area of 0.5 cm 2 , as is the case for the 


a cross 


Cary Model 82 Raman Spectrophotometer, then the required 


air flow rate is 24 liters^min 


-1 
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Suppose that we wish to use ambient outside air as our 
sample source, drawing it through the spectrophotometer by 
means of a vacuum pump as shown in Figure 9. We desire to 
minimize impaction of particles on the walls of the tubing 
leading to the spectrophotometer. We examine first the 
range of tubing diameters required if the flow through the 
tubing is to be laminar. 

The flow rate through the tube is related to the dia- 
meter of the tube and flow velocity via 

Q k7fd 2 V (55) 


We require that the Reynolds number for the flow be less 
than the critical value above which the laminar flow reg- 
ime gi'tes way to turbulence. According to Haltiner and 
Martin (1957), the Reynolds number is given by 


Re = 


v d 
max 


(56) 


Combining our two requirements, we find that 

v 

, 40 max 

Q SS 1 || I Cjj. » «m H- |» 

'TT^l Re v 


( 57 ) 


For air at a pressure of one atmosphere, in the atmospheric 
range of temperatures, 

H = 0.14 
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Sample Compartment 


Sample 

Cell 



V F 



Vacuum Pump 



Figure 9. Arrangement of pumping system for drawing 

particulate matter through spectrophotometer sample cell. 


Haltiner and Martin give the critical value of Re as 24,000. 
If we make the conservative assumption that 


then 


v 

max 

v 


crit 


- 7.6 x 10~Q 


If we take Q as 1000 co^sec** 1 , considerably greater than 
the minimum value determined above, we find that 

' d crit * °* 76 “ 

While considerations of wall effects limit the minimum 

• 

tubing diameter, the maximum diameter is limited by the 
sedimentation rdte of the particles; that is, if the velo- 
city of the particles through the tube is too small, they 
will sediment out of the air stream before passing through 
the entire length of the tube. 

If we assume that both the horizontal and vertical 
components of the particle motion are constant, we may 
write 

d - 2 2 L (53) 

v 

as the minimum allowable diameter for a given tubing 
length, for a particle entering along the center line of 
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the tube-. If we assume that the particles are spheres and 
that their terminal velocities are given by Stokes* Law 
(Fleagle and Businger, 1963), then 


w 


1 r ^ Part 
’ n f Air 


(59) 


The relationship, 

Q « k'Ttvd* 


(60) 


must be satisfied, as before. So, 

9Q * P Air 


If r Z p 


Apart 


gL 


(61) 


Assuming yalues for Q and >2 as before, and r ~ 10~^ CIU * 
^Air” 3 * 2 x 10 ' gra-cm" 3 , Ppart^ 2 and 


-2 

g - 980 cm~sec.* • we obtain 


d = 


4.91 x io 


For a 35 foot (1067 cm) length of ' tubing, the maximum 
permissible diameter is -thus 

d - 46 cm 

* 

Since, for a given flow rate, gas velocity increases as 
1/d and fall path increases as d, then tubing of diameter 
smaller than the value given above should be satisfactory 
from the standpoint of sedimentation. 
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* 


It thus appears that any tube size with 

0.76 cm < d C 46 cm 

is compatible with our flow requirements. The only re- 
maining consideration concerns the impaction of particles 
due to bends in the tubing. Some bends are necessary, but 
every attempt must be made to minimize their effect. 

(2) Execution of the experiment* 

The experimental arrangement employed here consisted of 
components essentially identical to those described above, 
including teflon tubing with a diameter of 1.9 cm. 

A schematic diagram of the area surrounding the 
spectrophotometer is shown in Figure 10. 

(a) Particle counter. 

In order to ascertain the extent to which the tubing 
affected the observed particle counts in the sample cell, 
a.Bausch A bomb* Model 40-1 Dust Counter was employed to 
measure particle concentrations at both ends of the tubing. 
The instrument is capable of counting only particles with 
radii greater than about 0.3M. It was felt, however, that 
if the tubing could efficiently transport particles of this 

wrom i - iwuxrfir ■ h i >,i .t, ,i ^ p.t«h ■jwiiiwtm ^v. j rn*-uuu'..,..-i nmi. ■' 

*Bausch & bomb, . Inc. , Rochester, New York. 




Figure 10. Scale diagram of area surrounding the 
spectrophotometer* showing the manner in which 
outside air was drawn in. 
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size, it would also transport smaller particles. 

The counter utilized a small pump to divert a small 
sample of the air flow through an optical particle counting 
apparatus. In attempting to use the counter, however, it 
was found that the pump was not sufficiently powerful to 
counteract the pressure drop produced when the air was 
pulled through the tubing. This, pressure drop was found to 
amount to about 35 mm Hg while, upon further investigation, 
it was found that the counter could not be operated at 
pressures differing by more than 0*22 mm Hg from the am- 
bient pressure of the instrument. 

In order to circumvent this requirement, the entire 
instrument was enclosed in an airtight piexiglas box which . 
was maintained at the pressure of the air inside the 
tubing (see Figure 11). 

Testing of the complete system, composed of tubing, 
sample cell, spectrophotometer, pump and dust counter, was 
Conducted at the Wright Chemistry Laboratory of Rutgers 
University, the Location of the spectrophotometer. At the 
time of the experiment, July 10, 1972, between 2:30 and 
3:30 PM ED ST, the counter indicated the following values: 






Figure 11. Diagram of enclosure that was constructed to maintain particle 

counter at pressure of air at the entrance to the spectrophotometer sample cell* 


■“4 
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Size 

Number Density (particles-cm **) 

Range 

Before Tiibing 

After Tubing 

Model used in 
Computations 

>0.3 m 

30 

30 

19 

>0.5^< 

4 

4 

4 

>1.0^ 

0.8 

0.3 

0.5 

> 2 ,0m 

0.0 

0.0 

0.04 


These observations are consistent with the model Intro- 
duced in Chapter IIX,H.3.b. In this and in another trial 
run conducted at New York University, it was apparent that, 
at least for particles with radii less than one micron, 
the tubing did not introduce severe losses. 

On the above* date, observing stations of the New Jersey 

Environmental Protection Department at Perth Amboy and at 

-3 

Somerville recorded particulate loadings of 45 Mg~m and 
~3 

123 respectively, values that are representative of 

urban areas (Junge, 1963)* 

(b) Sample cell. 

The Cary Model 82 Raman Spectrophotometer comes 
equipped with a roughly cylindrical multipass gas sample 
cell with Brewster angle windows at each end, and a single 
stopcock-equipped opening for filling. Since, for the pur- 
poses of this experiment, a continuously flowing sample was 
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required, a sample cell was purchased from Cary and modi- 
fied to add a second opening and to replace the existing 
opening with one that did not contain a constriction. 

Figure 12 shows, schematically, the original and modified 
gas cells. The cell, as modified, was found to introduce 
a pressure drop of approximately 115 mm Hg into the flow. 

(c) Vacuum pump, 

A pump was desired which could give a flow rate of at 
3 —I 

least 1000 cm -sec (2.1 cfm) at a pressure drop of 
150 mm Hg (35 mm in the tubing plus 115 nun in the sample 
cell). The pump selected was the Gast* 0822-V103-G271X, 
which is rated at 5.6 cfm under the above conditions. 

(d) Spectrophotometer, 

With all of the other equipment in place, experimenta- 
tion with the spectrophotometer controls indicated that the 

best possible spectra could be obtained with a spectral 

—1 —I — 1 

slit width of about 2 cm- and a scan rate of 0*4 cm -sec » 

These paramet£rs are somewhat less favorable than those used 

in demonstration air spectra shown in the manufacturer^ 

advertisements. The difference is attributable, first, to 


*Gast Mfg. Corp., P. 0. Box 97, Benton Harbor, Mich 49022. 



SAMPLE CELL 
BEFORE MODIFICATION 

Figure 12, Scale drawing of the original a: 

r • 


SAMPLE CELL 
AFTER MODIFICATION 


modified multipass gas sample cell 
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Che fact that the optical properties of the sample cell 
were, of necessity, damaged when the modifications were 
made. After modification, there was evidence of distor- 
tion of the Brewster angle windows which form the ends of 
the cell. In addition, the side walls of the cell, through 
which the Raman-scattered radiation passes, were slightly 
distorted. The second source of degradation of the Raman 
spectra was the fact that attachment iOf the tubing to the 
sample cell caused some slight dislocation of the cell, as 
well as backlash in the adjustments, for which no compen- 
sation could be made. 

(3Y Results * 

Several spectra were taken with aged, air conditioned, 
room air in the sample cell. One of these spectra is 
shown in Figure 13. The spectra were all identical, show- 
ing no peaks beyond those attributable to and 0^. 
Similarly, spectra taken with the sample cell removed en- 
tirely (un-aged, air conditioned room air containing sus- . 
pended particulate matter) were identical. 

Two spectra were taken with the pump in operation, 
dr awing, out side air through the sample cell (see Figure 
14). These spectra also contained no peaks beyond those 
found in clean air. 
























matter. X 
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A set of experiments was conducted to determine whether, 
perhaps, the return from particulate-laden air might contain 
a continuum of intensity across the Raman band that was not 
evident from observations of the individual line locations. 
Several pairs of observations were made, in each case with- 
out disturbing the adjustment of the spectrophotometer be-? 
tween observations. The first observation of each pair was 
of aged, air conditioned room air, while the second observe- 
tion was of outside air, with the vacuum pump operating as 
described earlier. The object of the experiment was to 
show that there is no correlation between weather conditions 
and the rotational Raman backs cat ter cross section of air. 

Since the indoor and outdoor temperatures did not dif- 
fer greatly, and since our attention is focussed on the 
least temperature-sensitive portion of the Raman spectrum, 
the assumption could be made that. the ratios between cor- 
responding 'individual lines in the two spectra were the 
same for any one case; that is, if the line intensities 
in the "outside air” spectrum were multiplied by a uniform 
factor, the resulting intensities would equal those oh? 
tained in the "room air" spectrum. To obtain the best 

estimate of this factor, about fourteen prominent peaks 

«*1 —1 

In the frequency interval, ll cm 98 cm were 
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Identified in each of the spectra, and the least squ^es 
regression line between the "outside air" and "room air" 
spectrum was computed. One such regression line is shown 
in Figure 15. The slope of this line is the factor that 
we are seeking, and its y-intercept is difference 

between the average background counts present 
spectra. The computed background count values were 
integrated over the entire spectral interval, 

11 cm -1 ^ Al> 4 98 cm" 1 and the resulting values normalized 
by; the total (signal plus background) counts obtained over 
that same interval. The resulting value is the "best fit" • 
differential background intensity over the interval, ex- 
pressed as a fraction of the total intensity. These values 
are shown in Table 5. The values shown are, in some cases, 
several percent of the total Intensity, and, in one case, 
almost twenty percent of the total. The computed differen- 
tial background intensities take on both positive and nega- 

tive values* 

No definite explanation could be found for these fluc- 
tuations. There is a possibility that small amounts of 
room light may have been ducted into the sample chamber by 

means of the translucent tubing used to bring the sample 

e i *r &ht could have both temporal 
air into the chamber* This lignt coui 



INTENSITY (OUTSIDE AIR) 




Table *>• Synopsis of measurements of meteorological parameters and concurrently observed Raman spectra. 
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1315 

301.1 

12.2 

. 71 *3 

2.0 

2t .1*5 

26.27 

0.727 

+2.358 

10.1*2 

♦0.199 

8/18/72 

1030 

29U.!i 

15.1 

72.0 

2.0 

1*5*56 r 

5U.5B 

0.838 

-0.162. 

10.01* 

- 0 . 01 U 

8/28/72 

1130 

298.3 

13.1 

75.6 

1.9 

7UJ*8 

95.82 

0.777 

♦0.012 

15.17 

♦ 0.0007 

8/29/72 

09US 

297*2 

13.0 

96.8 

1.8 

1(5.72 

61 .38 

0.753 

- 0 . 52 U 

9.59 

-0.0l*3 

9/13/72 

1000 

290.5 

13.6 

172. 

1.5 

1*5.36 

60.00 

0.71*6 

♦0.611* 

9*00 

♦0.060 

9/1 U/72 

1000 

297.? 

16.1* 

ioU. 

1.8 

113.85 

150.30 

0.756 

♦0.290 

8.70 

40.029 

9/15/72 

0900 

291 .6 

9*60 

113. 

1.8 

115.1*0 

1U8.1U 

0.777 

-0.1*29 

9.U6 

-0.01*0 

9/20/72 

0930 

287.2 

8.1*7 

75.9 

1.8 

HU.85 

11*6.68 

0.783 

-0.0$0 

8.88 

- 0.005 

9/22/72 

1030 

2S5.5 

11.8 

60.7 

1.8 

117.60 

153.12 

0.771 

♦0.369 

9.02 

♦ 0.036 

10/ 2/72 

lOUS 

291*1 

10.3 

89.9 

1.8 

113.55 

11*8.63 

0.765 

+ 0.357 

9.51 

+ 0.033 


* Prom hygro thermograph and high- volume sampler observations made at New Brunswick, N.J., about two miles from 
the location of the spectrophotometer . 

«* Signal levels and ratios indicated on the spectrophotometer were critically dependent on the slit width setting* 

' Ibis setting could not be. reproduced accurately enough to be considered as constant from one run to the next. 

Consequently 9 the slit width was set to the smallest value for which reasonably noise-free signals could be <x> 

obtained. ' ^ 
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and spectral variability. Furthermore 9 there is the pos- 
sibility of small fluctuations in the pumping rate. 

The thrust of this examination, however, was to deter- 
mine whether meteorological parameters influence the re- 
sults, At the timer of each set of observations, concurrent 
measurements of temperature, humidity and particulate mass 
loading were made, using a hygrothermograph and high volume 
sampler at the Rutgers University weather observation site, 
located about two miles from the spectrophotometer. A 
synopsis of the measurements is given in Table 5. Correla- 
tion coefficients were computed between the background 
values and the absolute humidity and the particulate mass 
loading. They are as follows: 

Differential background vs. Absolute humidity: +0.279 

Differential background vs. Particle mass loading: -0.157 
According to the significance test for correlation coeffi- 
cients given by Brooks and Carruthers (1953), neither of 
the. above Values is significantly different from zero. 
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Chapter IV. 

TEMPERATURE DEPENDENCE OF ROTATIONAL RAMAN LINE INTENSITIES 
AND DESIGN OF FAMILIES OF FILTERS WHICH NEUTRALIZE TEMPERATlff 
DEPENDENCE OF SCATTERING FROM THE ROTATIONAL RAMAN LINES OF.A 


List of Symbols 


- 1 . 


D 


C 

D 

d 

E(h) 

e 

e 

o 

G 

H 

h 

*P 

i: ' 


Doppler broadening (h&lf) half width (cm )% 

Pressure broadening (half) half width (cm ). 

Constant (cm~*-°K. ^). 

—1 o 

Constant (cm ~ K )* 

-1 

Speed of light (cm-sec ). 

Distance between interferometer and multipass 
reflection mirror (era)®. 

Width of a diffraction grating groove (era). 

Fraction of available Raman energy used by a filter 
with (half) half width h. 

Vapor pressure of sodium .in absorption cell (dyne—cm 

“2 

Reference vapor pressure of sodium (dyne-cm )» 
Gaussian filter form factor. 

Size of spot of light on interferometer (crn)» 

Filter (half) half width (cm *). 

Relative intensity of Raman line at frequency 
Number of temperature values used in determining 
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j 

K 


L 

M 

A 

m 

N 


n 


A 

n 



P 

p.nu 
r 
P„ 
s « HI© 


P 

R 

R* 

T 


mean line intensity. 

Rotational quantum number of initial state. 

Order number of principal energy maximum in light 

diffracted by diffraction grating. 

Latent heat of vaporization of sodium (erg-gm ). 

**1 

Atomic mass of sodium (gm~mol ). 

Refractive index of interference filter. 

Total number of grooves on diffraction grating, 
f /Number in an optical system. 

Number density of sodium molecules in absorption 
cell (cm ^). 

Number of passes of light beam through interferometer , 
Power incident on diffraction grating (W) . 

Power received in primary maximum of grating 

spectrum (W) . 

Value of P due to Raman scattering (W)* 
p.m. 

Power received in interval u, m £ n ^ ut ^ U *max of ■ 

secondary diffraction pattern of grating spectrum (W). 

Value of P due to elastic scattering (W). 
s.tn. 

Pressure (atm) • 

Reflectance of interferometer mirror. 

«“1 Q 

Universal gas constant (erg-mol - K ) . 

Temperature (°K)* 



89 


T Transmission of filter at center of its transmis- 

c 

sion curve. 

Reference temperature (°K) • 

Transmittance of absorption cell at wavelength X • 
Ty(h,lj) Transmission at frequency V, of filter centered at 
V C 3 with (half) half width h. 
tT Diffraction grating parameter relative to m=0. 

U 1 Diffraction grating parameter relative to arbitrary 

order number , m . 


U* 

max 


Kin 


W 


w 


Ai rms 


AL 


AV 



Value of grating parameter U' at 
Value of grating parameter U* at 9 mitl - 
Diffraction grating parameter. 

Width of interferometer surface (cm). 

Length of light path through absorption cell (cm). 
RMS fluctuation of 2Z about 2^ iy j expressed 
as a fraction of Xly « 

4V 

Geometric distance of diffracted Raman energy from 
diffracted elastic energy (cm). 

Free spectral range (&). 

Frequency interval over which Raman line intensities 
are summed (cm *). 

Deviation of caused by deviation of incident 

-1 

light from the perpendicular (cm ). 
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max 

^min 


e T 




Spectral distance of frequency for which 6=*9 from 

max 

that where 8~8 (cm" 1 ')* 

Distance between adjacent diffraction grating 
grooves (cm). 

Maximum allowable angle of light incident on filter- 
ing device. 

Diffraction angle of elastic scattering. 

Maximum diffraction angle of interval in which 
Ramaft lines are observed. 

Minimum diffraction angle of interval in which 
Raman lines are observed. 

Angle of light incident on diffraction grating, 
relative to normal. 

Angle of light diffracted by grating, relative to 
normal* 

o 

Wavelength (A) . 

* 

Frequency (cm ). 

Frequency at which filter transmission curve is 
centered (cm ) . 

Frequency of elastic scattering (cm ). 

Free spectral range (era" 1 ). 

Center of spectral interval V + $1/4^4 2 V - V ' ~ & V 

w V 

or 6 < ®min* 
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V r Frequency of Raman line, or nominal frequency of a 

—1 

group of Raman lines (cm ) . 

^ Object or image radius in an optical system (cm) . 

X.Iy Sum of all Raman line intensities at temperature T 

in frequency interval, AV, 

% I y Average of *5Lly over i temperature values* 

0^ Absorption cross section of sodium vapor at wave- 

2 

length X (cm ) . 

*7? Optical thickness of absorption cell walls* 
c 

(B Angle of tilt of interferometer relative to multi- 

pass reflection mirror. 
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A. General. 

As was mentioned earlier, although the total energy 
scattered in the 0- and S-branches is independent of tem- 
perature, individual line intensities do. vary considerably. 

In accordance with (20) and (21), as temperature is in- 
creased, the lines closest to the incident frequency de- 
crease in intensity, while those in the far wings increase. 
Table 6 shows some typical intensity variations to be ex- 
pected in the range of atmospheric temperatures. 

Obviously, depending upon which spectral interval of 
the Raman band the observer* chooses to monitor, rotational 
line intensities could be used as sensitive temperature 
indicators or as temperature-independent parameters . 

For the purposes of this study, it Is desirable to 
locate the family of spectral Intervals in which there is 
an insignificant net temperature dependence. Since the 
rotational Raman cross section, as defined in Chapter. III.F. 
is independent of temperature, this family includes, but is 
not l imi ted t6 the interval encompassing all lines (i.e., 
-<»<(y ~ V e ]< *») . However, we prefer not to use an eX 7 

tremely wide interval, since, as will be shown later in 
this chapter, background radiation would become a problem. 
Furthermore,, for practical reasons, we desire to examine 



Table. 6, Temperature dependence of some representative 
lines in the rotational Raman spectra of nitrogen and 


oxygen. 
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Av 

, "1\ 
(cm ) 

Midrange 
intens ity 
(re 1. units) 


0 

11.94 

9.276 

«2 

9 

83.51 

12.585 


23 

194.3 

0.256 


1 

14.39 

11.918 

°2 

11 

71.89 

20.335 


33 ' 

197.3 

0.049 


Percent change 
of intensity 
per +10°K 
temperature 


- n 
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intervals which lie only to one side of the elastic line* 

In the atmospheric range of temperatures, decreases in 
line intensity with increasing temperature occur only in 
the frequency intervals, 0 < jl/ r *v 70 cm*" 1 . Thus, 

any group of lines whose total intensity is to be indepen- 
dent of temperature must take in lines on both the high 
frequency side and the low frequency side of Jv r - V^j-70 cnT 


-1 


B. Ideal filters. 

For the purpose of illustrating the nature of these 
intervals, it was initially assumed that an ideal filter 
could be applied to the receiving apparatus; that is, it 

was assumed that one could choose to receive all the ener- 

* 

gy from any interval of lines within the 0- or S -branches, 
and exclude all energy from outside that interval. When 
this assumption is made, there can be only as many groups 
of lines with near-zero temperature dependence as there are 
lines in the region, 0 - V,| 70 cm 1 . There are, 

in fact, about eleven or twelve such lines in each branch. 
Intervals for near-zero temperature dependence were compu- 
ted, beginning with each of the first eleven lines in 
turn. The method of computation was as follows: 

Suppose that we are examining the total Raraan-scat- 
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tered intensity due to all of the lines lying in some 
frequency interval, AV , 

^ELly 

The average intensity observed over some temperature range, 
AT=T 2 -T i is given by 



Since ^>.1*, is a slowly-varying function of T, we may 
approximate the above integral by an average of ^Xv 
over i uniformly spaced temperature values in the interval 
of interest. That is. 



In our ideal filter we wish to minimize the temperature 
dependence of Sly • For the purpose of our computations , 
the interval of "minimum temperature dependence” was de- 
fined as the spectral interval, AV, in which the RMS 
deviation of JlI* from is a minimum, for tempera- 

ture variations over AT. That is, the interval, AV, for 


which 
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= minimum (62) 

i -i.lv 

Intensities were computed for each 20°K in the interval 
200°K ^320°K and are shown in Tables 7 and 8. The 

results of the calculations, shown in Table 9, indicate 
that, with an ideal filter, it is theoretically possible 
to achieve less than one percent temperature dependence 
while utilizing nearly all of the energy present in the 
0- and S -branches. 



C. Gaussian interference filters. 

1 . General . 

While maintaining the same criterion for "minimum 
temperature dependence" as was discussed above, we now 
extend our calculations by choosing a filter whose trans- 
mission function is Gaussian in character, with a spesi« 
fied center frequency and (half) half width (spectral 
distance from center frequency to the frequency at which 
transmission falls off by & factor of 0.5). That Is, 

(v- v \ 2 

„ ln2 \~H~-J 

G (h, V c ) = e (63) 
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Table 7. Computed line locations and relative intensities in 
the O-branch of the rotational Raman spectrum of air for 
temperatures in the range 200°K ^ T ^ 320°K. 
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7. si; 

1 . 656-03 

2.032-03 

2.44 

5- So 3-04 

7.163-04 

O.04 4 

2.066-04 

2.950-04 

4 , 0 :' 

T.5ftS-04 

9.086-04 

1*25.” 

9*C ll-OS 

1,370*04 

1.96 

2.457-04 

3.340-04 

4.SM 

3.693-6S 

5.983-05 

9,0£ 

3. 147-04 

4.420-04 


9.859-05 

' 1.432-04 

1.9? 

1 ,42**05 

2,466-05 

3.9? 

|.2p9-04 

1.819-04 

2.59 

5.154-05 

9.595-06 

1*04* 

3,6 ? 4-05 

5.6S8-03 

0.32 1 " 

4.255-05 

6,899-05 

1,04 

1 . 757-06 

3.S27-06 

6,4&: 

1,?E**G5 

2.061-05 

3.24: 

5*640-0 1 

1.225-06 

2,406 

1.376-C.5 

2.414~05 

3.92c 

3.783-66 

6.931-06 

1.17- 

1.7o»*07 

4,026-07 

MV 

4.077-r£6 

7.003-06 

1.371 

4,3s8*Cd 

1.251*07 

2.649 

lo O 70 -O 6 

2.153*06 

3.932 
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Table 8. Computed line locations and relative intensities in 
the S -branch of the rotational Raman spectrum of air for 
temperatures in the range 200°K 4 T ^ 320°K. 





Relative 





3 

Fre<tue^cy 

200°K 

220 °K 

0 

it 

0 

1*194+01 

1.141-c2 

1.038-02 

1 

0 

i 

1.439+01 

1.462-nZ 

1.332-02 

2 

N 

t 

1.990+01 

9.9R0-f>3 

9.100-03 

3 

0 

3 

2*690+01 

2.441-C2 

2.2 44-02 

<f 

M 

2 

2.705+Ul 

2. 693-02 

2.460-02 

5 

• H 

2 

3.501+01 

1.6G2-j,2 

1.479-02 

6 

0 

5 

3*740+01 

2.945-f,2 

2. 755-02 

7 

N 

4 

4.377+01 

3.506-02 

3. 272-02 

a 

0 

7 

4.891+01 

. 2.95a-(»2 

2.827-02 

9 

N 

5 

5.172+01 

1.409-i;2 

1 . 702-02 

10 

N 

6 

5.467+01 

3.503-02 

3.3h4-02 

u 

0 

9 

6.040+01. 

2.566-02 

2.539-02 

12 

N 

7 

6.762+01 

1 »626-(;2 

l.~?l-02 

13 

0 

11 

7.169+01 

1.979-»i2 

2.1*38-02 

14 

N 

B 

7.557+01 

2.493-02 

2.809-02 

15 

0 

13 

0. 336+01 

1.369-C2 

1.47Q-02 

16 

N 

9 

a. 351*01 

1.236-r2 

1.264-02 

17 


10 

9.145+01 

• 2.035-C2 

2.136-U2 

Id 

0 

IS 

9.483+01 

6.556-03 

9. /50-03 

19 

« 

11 

9.939+01 

Q.o?6-r3 

S.721“fl 3 

20 

0 

17 

1.063+02 

4.850-(«3 

5.881-03 

21 

N 

12 

1 .073+02 

1.230-02 

1.379-02 

22 

N 

13 

1 .’153 + 02 

4.504-r3 

5.283-03 

23 

0 

19 

1.177+02 

2.502*1)3 

3.252-03 

29 

N 

14 

1.232+02 

6«568-(i3 

7.850-03 

25 

0 

21 

1.291+02 

1.177-03 

1.653-03 

26 

N 

15 

1. 311+02 

2.277-r3 

2.829-93 

27 

U 

16 

1.390+02 

3.057-r 3 

3^60-03 

2 d 

0 

23 

1.405*02 

5«063~i'4 

7.734-04 

29 

N 

1 ? 

1.469+02 

9.942-o'f 

1.346-03 

50 

0 

25 

1.519+02 

1.992-n4 

3.336-04 

31 

N 

16 

1.546*02 

t.253-ft3 

1.770-03 

32 

N. 

19 

1.627+02 • 

3.828-r4 

5. /06-04 

33 

0 

27 

1.O35+02 

7.1<U-c5 

1,328-04 

39 

N 

ZO 

1-/06+02 

4.53o-<,4 

7.121-04 

33 

0 

29. 

1. 746+02 

2.374-n5’ 

4.886-05 

3b 

M 

Zl 

1.705*02 

1.303-C4 

2.160-04 

37 

0 

31 

1.060+02 

7.20D-C6 

1.661-05 

3fi 

N 

22 

1 .064+02 

1 ,'453-f 4 

2.550*04 

39 

N 

23 

1.943+02 

3.429- c $ 

7.310-05 

90 

0 

33 

1.973+02 

2.tl05-f,6 

5-285-Of, 

41 

N 

24 

Z. 021+02 

4.123-cS 

8.175-05 

42 

0 

35 

2.065*02 

5.132-07 

1.521-06 

43 

N 

25 

2 .100+02 

1.050-65 

2.222-05 

44 

N 

26 

2.179+02 

1.033-05 

2.350-05 

4b 

' 0 

37 

2.198+02 

1.203-07 

4,103-07 

46 

If 

27 

2.257+02 

2.493-{<6 

6.050*06 

47 

0 

39 

2.310+02 

2.614-oS 

1.025*07 

48 

N 

26 

2.335+02 

2.324-06 

6.064-06 

49 

N 

29 

2.414+02 

5.259-o7 

1.479-06 

50 

0 

41 

2.421+02 

5.206-09 

2.375-08 

hi 

N 

30 

2.492*02 

4 .623-97 

1.406-06 

52 

0 

43 

2.533+02 

9.546-lfl 

S.105-G9 

33 

N 

31 

2.570+02 

9.870*08 

3,251-07 


Relative 

Intensities 




gU0°K 

260°K 

280°K 


2* 

9 » 515-03 

8.786-03 

8.1^0*03 

7.617-03 

. 7.14: 

1.223-02 

1.130-02 

l«05l“02 

9.020-03 

9.21 

8.362*03 

7.735-03 

7.1o5-03 

6.726-03 

6.31* 

2.077-02 

1,933-02 

l.Co8-02 

1.697-02 

1.60- 

2.278-02 

2.115-02 

1.973-02 

1.050-02 

1.74 

1.374-02 

1,203-02 

1.2 0 3-02 

1.132-02 

1.06* 

2.586-02 

2,435-02 

2,300-02 

2.179-02 

2.07' 

3.066-02 

2.884-02 

2«7?i" 02 

2.575-02 

2.4*) 

2.709-02 

2,595-02 

2-4n8-02 

2.308-02 

2.2.9; 

1.612-02 

1,530-02 

1.455-02 

1.306-02 

1*32. 

3.229-02 

3,098-02 

2.974*02 

2.857-02 

2.74". 

2.499-02 

Z, 449-02 

2.3g3*02 

2,335*02 

2.27/ 

1.550-02 

1.506-02 

1.462-02 

1.418-02 

1.37 

2.072-02 

2 . 088—02 

2.O9O-02 

2.002-02 

2*06. 

2 » R 65*?02 

2,025-02 

2.777-02 

2.723-02 

2,66'' 

1.563-02 

1.620-02 

1.677-02 

1.712-02 

1.73^ 

1.270-02 

1,281-02 

1.277-02 

1.268-02 

" 1.25' 

2. 206-02 

2.253-02 

2.2^2*02 

2.296-02 

• 2*29* 

1.079-02 

1.160-02 

1.243-02 

1.306-02 

1.35 k 

9.227-03 

9.616-63 

9.9o«*03 

1.012-02 

1*025 

6.853-03 

7.750-03 

8.5&5-03 

9-295-03 

9.944 

1.497-02 

1,595-02 

1.675*02 

1.739-02 

1-79'. 

5.900-03 

6,438*03 

6.6n9-03 

7.290-03 

7.614 

4.016-03 

4.770-03 

5.4 9 7-03 

6.187-03 

6. 83: 

9,037-03 

1,012-02 

l.l0S“02 

1.194-02 

1.269 

2.176*03 

2.729-D3 

3.205*03 

3*061-03 

4.41/ 

3.365-03 

3.072-03 

4.343-03 

4.774-03 

5.16b 

4.876*03 

5.778-03 

6.645*03 

7.466-03 

8.23T 

1.093-03 

1.454-03 

J. 84 0*03 

2.263*03 

2.6?r 

1.719-03 

2.102-03 

2.4p3-03 

2.055-03 

3.21; 

5.0H8-04 

7.226-04 

9»7(j6-04 

. 1.248-03 

1*540 

2.361*03 

2,983-03 

3.625-03 

4.271-03 

4 0 9 1 1 j 

7.897-04 

1.033-03 

l.2 9 3-03 

1.564-03 

1.03C 

2^201-04 • 

3.352-04 

4. 7^2— 04 

6.474-04 

0.40^ 

1.029-03 

1,396-03 

1.B04-03 

2.241-03 

2.69'T 

fl. 048-05 

1.453-04 

2.2,1 1-04 

3*166-04 

4.316 

3*267-04 

4.607-04 

6.150-04 

7.863-04 

9.70L 

3.309*05 

5.090-05 

9.6 0 3-05 

1.460*04 

2.097 

4.044-04 

5.936-04 

8.2„a-04 

1.081-03 

1.370 

1*220*04 

1.860-04 

2.676-04 

3.637-04 

4,73-7 

1.152-05 

2.234-05 

3.920-05 

6.351,-05 

9.640 

1.435-04 • 

2.296*04 

3.4i6-04 

4.796*04 

. 6.42f 

3-735-06 

7.933*06 

l.So5-05 

2.609-05 

4.204 

4.1 17-0S 

6.694-05 

1.066-04 

1.549-04- 

•—2.139 

4,607-05 

0.089-05 

1.3r4-04 

1.962-04 

2.793 

1.128-OS 

2.639-06 

5. 4j6- 06 

1.012*05 

1.736 

1.257*05 

2.319*05 

3.800-05 

6.085-05 

a,94e 

3. 178-07 

8.226-07 

1.848-06 

3.711-06 

6.799 

1.338-05 

2.599-Q5 

4*564-05 

7.401-05 

1.125 

3.476*06 

7,117-06 

i.3o8-05 

2.206-05 

3.471- 

8*351-08 

2.404-07 

5.919-07 

1.286-06 

2.525" 

3.424-06 

7,620-06 

1.4&8-05 

2*578-05 

4.203- 

2.047-08 

6,591-00 

1.706-07 

4,215-07 

8.900- 

9.714-07 

1.994— Q6 

4.0J2-06 

7 e 386-06 

1.249- 



Table 9 ft Characteristics of various ideal filters. Line numbers refer to 
consecutive numbers assigned to each line in Tables 7 and 8. 


First Line 

LLast Line 


(Half) 


Percent 

Line 

Line 

Line 

Line 

Center 

Half 

Fraction of 

RMS Intensity 

No. 

Freq. 

(cm- 1 ) 

No, 

Freq. 

( cm"’-*' ) 

• Freq. 
(cm~ s ) 

Width 

(cm“^) 

Available 
Energy Used 

Variation over 
200°K~320°K 





0 - Branch 



1 

14.39 

24 

123.2 

68.79 

54,4 

0.9483 

.0.38 

2 

19.90 

22 

115.3 

67.60 

47.7 

0.8609 

0.54 

3 

25,90 

21 

107.3 

66.60 

40.7 

0.8276 

0.59 

4-.- 

27.85 

20 

106.3 

67.08 

39.2 

0.7453 

0.69 

5 

35.81 

18 

94.83 

65.32 

29.5 

0.6551 

0.82 

6 

37.40 

17 

91.45 

64.43 

27.0 

0.5995 

0.85 

7 

43.77 

15 

83.36 

63.57 

19,8 

0.4676 

0.91 

8 

48.91 

14 

75.57 

62.24 • 

13.3 

0.3622 

0.98 

9 

51.72 * 

13 

71.89 

61.81 

10.1 

0.2403 

0.99 

10 

59.67 

12 

67.62 

63.65: 

3.98 

0.1594 

1.20 

U 

60.40 

11 

60.40 

60.40 

0.0 

0.0555 

mm 





S - Branch 



4- 

27.85 

33 

163.3 

95.5$ 

81.7 

0,8979 

0.23 

5 

35.81 

28 

140.5 

* 88.16 

52.4 

0.8305 

0.40 

6 

37.40 

26 

131.1 

84.30 

46.9 

0.7833 . 

0.51 

7 

43.77 

24 

123.2 

83.49 

39.7 

0,7107 

0.62 

8 

48.91 

22 

115.3 

82.U 

33.2 

0.6086 

0.78 

9 

51.72 

20 

106.3 

79,01 

27.3 

0.4972 

0.79 

10 

59.67 

19 

99.39 

79,53 

19,9 

0.4435 

0,89 


vO 

y0 
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Obviously, there will be an infinite number of combi- 
nations of center frequency and half width that will satis- 
fy our requirement for a small temperature dependence. We 
recognize immediately, from 4n examination of Table 9, that 
for a filter of near zero half width, the desired center 
frequency will fall at about 60 cm*’ 1 in the O-branch and 
75 cm** 1 in the S-branch. Furthermore , we know that as the 
half width is increased, the desired center frequency will . 
move progressively outward into the wings of the 0— and S— 
branches. 

Calculations were made of the various parameters of 
interest for filters with half widths of 25, 37.5, 50 and 
75 cm”* . Plots of net temperature dependence for various 
center-frequency /half width combinations are shown in * '• 

■ ’ * v. ^ 

Figure 16, illustrating the characteristics discussed above. 

From the results shown in Figure 16, it is a simple 
matter to interpolate the center frequency /half width 
combinations that result in minimum intensity variation. 
These combinations, together with limits for <1% tempera** 


ture dependence, are shown in Figure 17. The fraction of 

available Raman energy used, E(h), is given by 


E(h) 



(64) 


0 



Percent RMS Temperature Variation Over 200 9 K ^T^320® 



Center Frequency of Filter Relative to Transmitted Frequency (cm- 1 ) 

Figure 16, Net temperature dependence of energy transmitted through a Gaussian. * 

filter for various center frequency /Kalf width combinations. 

. « *■ * . . , 
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Center Frequency and Limits for < {% intensity Variation (cm* 1 ) 



10 20 30 40 50 



(Half ) Half Width (cm- 1 ) 


Figure 17* Temperature dependence field for energy transmits 


10 


Here, h and V are ordered pairs from Figure 17 which give 

V- 

minimum temperature dependence. A. graph of E(h) versus h 
is shown in Figure 18. 

Since we find that there are an infinite number of 
center frequency /half width combinations from which to 
choose, we shall investigate, in the following sections, 
some additional constraints which may be applied to the 
filter design. . In the next section, we shall investigate 
the possibility of using the filter to suppress the signal 
due to elastic scattering from the Raman channel. In 
Chapter IV. A.', we shall consider the possible selection of 
filter parameters to optimize the signal to noise ratio. 


2. Collimation requirements. 

When using interference filters, or any interferomet- 
ric device for that matter, it is necessary to examine the 
requirements that the device imposes upon the collimation 
of the incident light. From the optical geometry of inter- 
ference filters, it can be shown that the incidence angle, 

4 

0, and center frequency deviation, Av> q , of a filter are 


related according to (Oriel Optics Corp. , 1971) 


sin © - m 



AVol 

V, 


h 


(65) 



Relative Filter Efficiency For 


1.0 
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0 10 20 30 40 50 60 70 80 


(Half) Half Width (cm* 1 ) 

Figure 18 * Fraction of available Raman energy transmitted by 
Gaussian filters of optimum center frequency and half width* 
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Typically, m, the refractive index of the filter, is about 
1.7. For the types of filters discussed here, we shall 
require the center frequency to vary less than about 5 cm 
with changing incidence angle. Then 

e < 2.5° 

Figure 19 shows the interrelationship between f , and 9 

in an optical system. In a system which introduces no 
losses or aberrations, it can be shown (Born and Wolf, 1959) 
that the quantity f iJf is constant; that is, the optical 
properties of the system may be transformed along lines of 
constant f / J/* . This is a convenient parameter for ex- 
pressing our results, in that it allows the comparison of 
systems which require different blur radii and f /numbers. 

A filter of radius 2 cm, with a maximum acceptance 
angle of 2.5° falls at the location marked (£) in Figure 19. 
We see that, for an optical system to produce such a set of 
parameters, it is required that 

fuM 3 < o - 17 

A further discussion of the practical attainability 
of such a system will be presented in the following sec- 


tions . 
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D» Interference filters which provide sufficient suppres - 
sion of the elastic scattering in the Raman channel, 

1 , General , 

It was shown in the previous section that, if in 
designing a filter system, we consider only neutralization 
of temperature dependence of the received signal, we have 
an infinite number of filter transmission functions from 
which to choose. With this fact in mind, we examine the 
possibility of additionally taking advantage of the fil- 
ter 1 s . transmission function to suppress the elastically 
scattered light from the Raman channel. 

In Chapter III.G; it was shown that S-branch Raman 
scattering is 



-30 

2.85 x 10- u 


2.01 x 10 




- 0.0142 


times as intense as the elastic molecular scattering, and 
that the 0-branch Raman scattering is 0.0102 times as 
intense o . • . 

We desire that the elastic scattering be suppressed 

to less than one percent of the Raman scattering, in the 

* * * . 

Raman channel. This implies a required transmission at the 

-»4 

elastic frequency 1.42 x 10 times the effective Raman 

-4 

band transmission, E(h)T_, and a .factor of 1.02 x 10 for 
* c 
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the 0~branch. When the effect of particulate elastic scat-: 
tering is included, we find that the required rejection 
factor not only becomes larger, but also wavelength^depend- 
ent* In the spectral interval 0.35 JU 7 0.70^, the inten- 
sity of molecular scattering varies by a factor of sixteen, 
while that of particulate scattering varies by a factor of 
less than two (Elterman, 1968). Thus, at 0.70X^s the re^ 
jection requirements for the filter are about ten times as 
stringent as at Q*35M* 

Elterman* s tables, which give attenuation coefficients 
as a function of wavelength for an atmosphere with a visi- 
bility of approximately 14 miles, indicate that the ratio 
of particulate to molecular attenuation coefficients near 
the surface is about 30 at and 3 at 0.35^ * 

Deirmend jian* s (1969) study of particulate phase functions 
indicates that the ratio of particulate to molecular back- 
scatter cross sections should be of the same order of mag- 
nitude as that of the attenuation coefficients* 

Under extremely hazy conditions, we might find these 
factors to be as much as five times as large (visibility 
less than three miles); that is, the ratio between parti- 
culate and molecular backscatter could be as high as 15 at 
0 * 35 ^/, and 150 at 0.70^ under these extreme conditions 
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and the rejection characteristics of the interference fil™ 
ter must accordingly be specified 15-150 times as high as 
for a purely molecular atmosphere* The particulate contri- 
bution to the backscattering diminishes rapidly with height 
the ratios of particulate to molecular backscattering are 
less by a factor often at 3 km, as compared to the ratios 
near the surface. Thus, for systems whose operation is 
confined to higher altitudes, the rejection requirements 
are not as stringent® 

Thus, we see that, for an interference filter to pro- 
vide the required suppression, the elastic frequency must 
fall far in the "wings’ 1 of the transmission function; 
that is, we require a relatively narrow filter whose trans- 
mission function drops steeply in the wings. In order to 
obtain design parameters for such a filter. Infrared In- 
dustries, Inc.*, a concern which manufactures interference 
filters was contacted. They supplied a set of typical 
transmission curves for filters that they produce. The 

filters are classified as one- through four-period filters, 

« 

the term "period" referring to the number of sets of di- 
electric layers composing the filter. Four pe.riods is 



*F. 0. Box 557, Waltham, Mass. 02154. 

* .*- 4> 


w* 


1X0 


presently the maximum number commercially available* 
Through use of the curves, which are not Gaussian in 
character (See Figure 20), it was possible to find the 
widest filter of this type that would give a rejection of 
10“ 4 or 10*" 5 at 6943 A, and still give the desired non- 
temperature dependent Raman transmission. This was done, 
using the same computer program as was used in the calcu- 
lations for ideal and Gaussian filters. The computed 
characteristics of such filters are shown in Figures 21 
and 22, and in Table 10* 


Table 10. Widest four period filters which give a rejec- 
tion of 10~ 4 and 10~ 5 at 6943 A, and still give the 
desired non-temperature dependent Raman transmission. 



Center 
tfave length 
(A) 

(Half) Half 
Width 
(£> 

Fraction 
of Raman Branch 
Transmitted 

, 

a 

o 

A U 

s o 
O 

'' 1 
& 

O^Branch 
S -Branch 

6212.9 

6980,6 

8.1 

10.0 

0.361 

0.384 

d 

O 

*r4 

st 

0-Branch 

6912.3 

11.7 

0.500 

I o 
O 

tH 

G> 

p4 

S-Brahcb 

6981.3 

14.6 

0.523 










.001 


1“ Period 
\ 

N» 


\ \ Z - Period 


W 


Gaussian 


3 -Period 


4- Period 
(estimated) 


.0001 


NUMBER OF (HALF) HALF WIDTHS FROM CENTER FREQUENCY 

Figure 20, Transmission curves. of typical single-* and multi- 
period interference filters (after Infrared Industries data) 




(HALF) HALF WIDTH ( cm-* 8 ) 


Figure 2£« Temperature dependence field for energy transmitted through 
four period filters of optimum center frequency and half width* 



iO 

(HALF) HALF WIDTH (cm- 1 ) 

Figure 22. . Fraction of available Raman energy used by four-per 
filters of ODtiraum center frequency and half width. 
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2. Collimation requirements* 

Obviously, filters of this sort require a high degree 
of collimation. As a first approximation we can assume, 
with theoretical justification, that the entire transmis- 
sion pattern of a filter is shifted in wavenumber space, 
without distortion, when the rays are incident at some small 
off-axis angle. Suppose we require that the off^axis rays 
may not change the center frequency (and, hence, that of 
the entire transmission pattern) by more than 2 cm- . For 
a filter centered at any of the wavelengths in Table 10, 
this amounts to a shift of less than 0,125 (half) half 
widths <. 

Such a shift has two effects that are of concern to us 0 
First, there is an alteration in the degree to which the 
filter balances the temperature-induced fluctuations in the 
intensity of the Rarrian signal, and second, there is ari 
alteration in the degree to which the elastically-scattered 
signal is rejected* 

The magni tude of the first effect may be determined 
through examination 6f Figure 17, and is found to fall 
within the one percent limitation lines. The variation in 
the rejection of the elastically-scattered signal may.be 
ascertained by inspection of Figure 20, At spectral dis- 
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tances beyond one (half) half width from the center fre- 
quency, the filter transmission variation induced by a 
frequency shift of 2 cm is less than a factor of two*. 

The off-axis angle corresponding to a frequency shift 
of 2 era is, according to (65)^ 

© < 1 , 6 ° 

Proceeding as in Chapter IV. C. 2., we find that, for a fil^ 
ter of radius 2 cm, it- is required that 

f tjf* < 0.11 

which falls at the location marked (b) in Figure 19. 

E» Filtering through use of a diffract ion_gr ating, 

1, General,. 

The possibility of using a diffraction grating to 
separate spatially the elastic and Raman components of the 

received signal was considered* 

While the necessary spectral and spatial.. resolutions 
are readily obtainable in commercially available gratings, 
the secondary diffraction pattern generated by the strong 
elastic signal is observed at the same angular position as 
the primary Raman lines, and can interfere. An investiga- 
tion is made here into the magnitude of this problem, and I 
into the design of a grating system which would minimize 
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it* 

Consider a fully-illuminated grating with N grooves. 
The energy diffracted per unit angle is given by Ditchburn 
(1963) as 


dF^ sin 2 U sin^NW 

de’ 1 U 2 sli?W~ 


( 66 ) 


where 


rrvd 


u » 


(sin9* * sinQ^) 


(67) 


and 


s/ £ 

W = « — (sin9* - sine) 

c 


( 68 ) 


If d = <£ ? as is approximately the case in modern gratings, 
then U » W, and 

..dP 9 sin^Ntf 

™ — j- . (69) 

de* 1 (NU) •• 

If the grating is blazed (Jenkins and White, 1957) such 
that the optimum order at is K s then (69) becomes 


dP. - stn 2 NU' 

— - - P.JT — , (U> s U - K'Tf) (70) 

d@> 1 (NU') 
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We define the "principal maximum" (p*m*) as the energy ly- 
ing between the first energy minima to either side of 
U — K'T?' (U 1 ~ 0), the center of the pattern* These minima 
occur at U* = * if/N (see Figure 23), and thus the energy 
in the principal maximum is 



r\H 

N dP d© 1 

— - — - du* 

• dQ* dU f 
if 

l/"S 


ss. 


2 cP.N 

V d cos© 1 
e 


if 

O 


sjLn 2 NU> 

* d(NU») 

(NU 1 ) 


(71) 


and the energy in some interval, ^ U* ^ U^ ay , o£ 

the secondary diffraction pattern is 


cP t N 


sin 2 ™ 


s*m* 


H'Vd cose* 

e 


IU* 
min 


(NU S ) 


d(NU* ) 


(72) 


2 » Constraints on the order number and free range of the 
gratin g sy stem* 

■■ ■ ■ urn *■ — i ■ii rffW a n ■ i l h m 

We require that the free range of the grating system 
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H 1 ■ ■ - I" -- ■■■.'■ —■I 

nS sin@ 0 , • sin 0 max sin 9 ml 


Figure 23.. Patterns due to elastic line and two Raman lines 
as generated by a diffraction grating (not to scale) . 
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be some minimum value, V^, so as ** nsure that the pat- 
terns due to adjacent orders do not overlap to an extent 
that would interfere with observations* Consider an ob— 

• servation made at an angle 0 Q * Using (68), we see that the 
principal maximum of frequency V Q diffracted in order K 
will be observed at this angle if 

sin 0 - sin 0 « (73) 

° •<* * e 

Similarly, the principal maximum of frequency V e + 
fracted in order K+l will be observed at this same angle if 


sin 9. 

. e 


sin 9 . 
o 


c (K+l) 

d ( y e + y f ) 


(74) 


Combining (73) and (74), we see that the criterion for both 
frequencies to be observed at the same angular location is 


that 



(75) 


This constitutes the maximum usable order number if £re~ 
quency V p is not to interfere with V e + The magni- 

tude of is determined by the character of the received 
signal and by the additional filtering in the system. For 
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instance, if we were certain that the received signal, in- 
cluding background, extended no farther than 1000 cm - * to 

either side of V , taken as 14,400 era (6943 A), then the 

6 

maximum order number would be 14,400/1000 ~ 14. Thus, a 
system of this sort must operate in the range 

1 < K ^ 14 

approximate ly * 


3. Interference of diffraction due to elastic scattering 
with that due to Raman scattering. 

Suppose that we are interested in observing the princi- 
pal maxima due to a group of Raman lines lying in the in- 
terval ( V t 4 V ’C (2 V ~ V - &V ) (see Figure 23). 
Then, applying the grating equation at these frequencies, 
we obtain 


sin & . - sin © 

min o 


cK 


d ( y + iv) 


(76) 


and 


sin 0 


max 


sin 0 = 

o 


d (2y - y • 

o e 


(77) 


Therefore, using (67), we find 


V'. ** Ktf . 

wtn 1 V„+ 


(78) 
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and 


U f » 
max 


K 


( zy 2 v^ \ 

\ 2 V 


(79) 


If the Raman energy received is presumed to vary with 

l/ + &y through E(h), as given in Figure 18, then the ratio 
6 

of primary Raman scattering to secondary elastic scatter- 
ing in the interval is 


^s*m. 

e 

P 

p.nw 

r 



(80) 


The value of the integral in the denominator is a constant, 
approximately equal to 4 lf/2. 

To estimate the values taken on by the integral in the 
numbfator for a given set of limits, we let 


stn^U 1 - sin 2 NU* = ?2~ /Z 
and. integrate; 
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. (l v -2 y - £v\ 

\ ° e d(NU') 


'Wkn <fy 


<V*i> -fr~ 


2 (y o -y) - 


(81) 


Thus, the ratio becomes, approximately. 


P 

V 


P 

P*BU 


r 


P e W 

P E(h) KN 
r 


< V q- V e) ~ ^ 

2< v o - y e ) 


(82) 


Where the wavenumber resolution is equal to the quotient, 
V /KN (Born and Wolf, 1959). The term .. 

w 


yr 

iy 2 Jv 


( v 0 - v e ) 

2 (v 0 -)2 e ) 


is determined by our previously described requirements, is 

independent of our choice for and takes on a value of 

approximately lcf 2 . The magnitude of the factor P e /P r E(h) 

is determined by the ratio of the elastic and Raman back- 

scatter cross sections and the fraction of the total Raman 

energy utilized, is also relatively independent of the sys- 

2 

tem parameters and is approximately equal to 1.7 x 10 . 

The ratio, P 0 m /P^ m * is determined by the level of ac- 
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curacy that we wish to achieve in our observation and, for 
one percent accuracy, is equal to 10 , Thus determined, 

the required resolution of the grating is found to be ap- 
proximately 0,006 cm , a value essentially independent of 
the system parameters. Such a grating, particularly one of 
good quality, is exceedingly expensive. It should be re- 
membered that, even though a grating is used to separate 
the components, filters are still necessary to produce the . 
desired weighting function. 

The required parameters for the grating may be relaxed 
somewhat by using a combination of grating and filter to 
produce the rejection. For instance, if a filter were em- 
ployed which had a transmission of 10 at then the 

resolution .of the grating system could be reduced by a 
factor of ten. 

4, Interference due to grating -"ghosts* 1 . 

The use of a grating presents another problem, one 
whose magnitude is not as readily determined. That is the 
problem of grating "ghosts 11 or anomalous diffraction 
properties which appear due to errors in ruling the grating. 
These patterns also tend to appear in the region where the 
Kaman lines are observed, and vary in location and inten- 


sity from one grating to another. 

5. Coll imat ion requirements. 

In order that we may focus the received rays closely 
enough to allow them to pass through the monochrome ter 
entrance slit, we must impose a constraint upon the blur 
radius produced by the optical system. 

The maximum allowable size of the entrance slit (and 
thus, that of the blur radius) is determined by the spec** 
tral resolution that we wish to achieve. in the monochro- 
meter exit plane, as well as by the optics comprising the. 
instrument. 

Fundamentally, we require that the monochrome ter be 
capable of separating the primary diffraction due to elastic 
scattering from that due to Raman scattering. That is, the 
dispersion of the instrument must be such that, in the exit 
plane, the distance, A L, of the Raman- scattered energy 
from the elastically scattered energy must be at least as 

great as the half-width of the entrance slit* The expres- 

* 

sion for the linear dispersion of a monochrometer may be 
written as 


KM J/° 



125 


If we parameterize this expression In terms of the blur 
radius of the optical system, F , and solve for the ratio 
F / Jf* discussed in Chapter IV 8 C*2., we find that 

4 - — <S4) 

of Al* V V cos 0 

vS w 

For an estimate of F we suppose first that the 
distance, in the exit plane, between the Raman-scattered 
and elastically scattered energy, Ah, is equal to the 
radius of the focal region, ]?• That is, that 

■**1 

FUrther, we assume that the resolution, 1^/KN, is 0*008 caT , 

compatible with the computations of the previous section* 

The frequency difference between the elastic and Raman 

—1 

lines, AV', is about 25 cm and cos 0 is not much dif- 

”1 

ferent from one* If V « 14400 cm , the ruby laser 

w 

frequency, we find that 


X 

•r 


0*22 cm 


This figure is comparable to those obtained when one 
mines the constraints on interference filters (see Figure 19} * 
That is, the constraints on the optical system utilizing a 
diffraction grating are of the same magnitude as those on . 
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the system employing an interference filter. 

F. Filtering through use of a Fabry-Perot interferometer. 

1. General. 

According to Jenkins and White (1957), the reflectance 
of a Fabry-Perot interferometer varies with wavelength in 
the manner shown in Figure 24. Regardless of the magnitude 
of R, the reflectance of the individual interferometer ele- 
ment, the reflectance of the interferometer always reaches 
zero at points spaced periodically in wavelength space. 

This is not true in practice, since the reflectance minimum 

• 

is limited by random scattering from the interferometer 
surfaces. The magnitude of this scattering is of the order 
of one percent o 

2. Constraints on the design of the inter ferometer* 

We wish to design an interferometer which has reflec- 
tion characteristics such that the region of near-zero 

reflection is wide enough to accomodate a laser line with 

o 

a spectral width of about 1 A, yet sufficiently narrow so 

o o 

that the reflectance in the region about 10 A. - 70 A may 
from the line is fairly large * 

o 

These requirements lead to selection of X 80 A, 
and mirror reflectances of about 25 percent* The thebre— 





128 


tical reflectance curve for such an interferometer, super- 
imposed on the rotational Raman spectrum of air, is shown 
in Figure 25. 

Since the free range is large, an interferometer of this 
design is not as sensitive to the angle of incident light as 
is a sharp transmission filter. Suppose that the mean in- 
cidence angle is five degrees (The angle must be as small 
as possible, in order to minimize polarization effects). 

Then, a one degree deviation in this angle will produce a 

o 

wavelength shift of 12 A in the interferometer pattern. 

This is much too large, of course. The sensitivity of an * 
interferometer to angular deviation of the incident light 
decreases as the angle of the incident light decreases. 

The limit is reached when the mean incidence angle is zero 
degrees. Suppose that this is the case, and that the beam 
has a one degree deviation to either side. In that case. 



cos 0^ 


A 2 

cos ©2 


A 1 

cos 

1° 

6943 A 

cos 



0.99985 


(85) 


o o 

Therefore, =* 6941.95 A or 13 1.05 A 








130 


Since this is about the maximum spectral deviation that we 
can tolerate, it is obvious that our system must operate 
at reflection angles very close to zero degrees* Suppose 
that the beam is centered at one degree to the normal* 

Then, 

^1 cos h° 

— — i i— m m m ss "**** ,B ■ 

6943 A cos 1%° 

«* 6945*08 A or A 2 ~ \ ^ B °S A 

This is about the minimum deviation that we can realistical- 
ly expect from a beam divergence of one degree* If we de- 
crease tlie divergence to 0*5°, then 

cos 

6943 A cos 1° 

~ 6943*76 A or ” 0.76 A 

So, if the beam divergence of the incoming light is 

10~ 4 radians, then the maximum ratio between the diameter 

of the collecting mirror and the Fabry-Perot interferometer 

is ^ 

0.5® / 57.3 



87*3 
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For an effective Interferometer spot 1 cm wide, the maximum 
collector diameter is 87.3 cm (34.4"). For a beam width of 
1C T 3 radians, this diameter is only 8.73 cm (3.44"). 

Suppose that we wish to improve the rejection charac- 
teristics of our system by allowing the beam of light to 
pass n times through the interferometer. The geometric 
aspects of such a system are shown in Figure 26. 

The effect of passing the beam through the interfero- 
meter n times is to produce a new overall transmission 
function that is determined by raising the original func- 
tion to. the nth power. This produces a curve whose rejec- 
tion characteristics are greatly enhanced, but whose 

transmission- qualities are mediocre. For example, an . 

-2 ‘ . 

interferometer : whose null reflectance Is 10 and whose 
"band pass" reflectance is 0.6, when used in a 3-pass con- 
figuration, will theoretically yield a null reflectance of 
iq~ 6 and a "band pass" reflectance of 0*22. 

Referring to the figure, we see that, in order to just 
separate the incoming and reflected beams, the distance be- 
tween the interferometer and the reflecting mirror must be 





where H is the size of the individual spot on the inter- 
ferometer, and that the minimum width of the interferometer 

must be 


w- 


1 + 


1 





(87) 


If e » fc°, Q * 3/4° and H ~ 0.29 cm, then 


•> . D •=» 17 cm. 

If n » 3, then 

w - 1*45 cm» 


3, Collimation requirements^ 

The value of f /Jf associated with an interferometer 

of this sort is about 0.004 cm, a value much smaller than 
those computed for the interference filters or diffraction 
grating. This is attributable to the fact that, in a 
multipass arrangement, it is necessary to keep the spot 
size verjr- small in order that the total size of the inter- 
ferometer be kept reasonable, and because the individual 
spots may not overlap. The interference filters discussed 
earlier are really nothing more than a large number of 
Fabry-Perot interferometers in a series or cascade arrange- 
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merit. The sharp pattern Is not produced by multipassing, 
and thus a much larger spot size may be used to obtain 
given interferometer pattern for a given interferometer 

size • 

It is to be concluded that use of a multipass Fabry- 
Perot interferometer results in severe restrictions on the 
optical quality of the total system, and that from this 
point of view, the other devices discussed earlier are 

highly preferable, 

c Filter system utilizi ng selective absoretionof 
elastical ly scattered light. 

1. General. 

In view of the high optical quality required in systems 
which attempt to separate, interferometrically, the strong 
elastically scattered light from the much weaker Raman 
scattering, an investigation was made of the possibility of 
selecting a laser wavelength which corresponds to an Iso- 
lated absorption line in some substance. A measured amount 
of that substance could then be introduced along the light 

path in the receiver, thereby absorbing the bulk of the 

elastically scattered light, and reducing its intensity 
to a low enough level that rejection would no longer be a 

problem* 
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2. Choice of an absorbing substance . 

In choosing an absorbing substance, the following 
factors were considered: 

(1) The absorption line should be isolated, with the 
nearest significant line more than about 150 cm 1 away. 

(2) The absorption line must be so narrow that it does 
not cause. significant absorption at the Raman frequencies. 

(3) The absorption line must be strong enough that 
transmission at its center can be reduced to \ 10 times 
that in the wings, using a manageable amount of absorbing 
material* 

(4) The transmission properties of the absorbing 
material must be such that defocussing of light at the 
Raman frequencies (e.g., due to spatial fluctuations in 
refractive index) be minimized. 

(5) The chemical properties of the substance must be 
such that it is not unduly dangerous or awkward to handle 
in the laboratory or in the field. . 

(6) The wavelength of the absorption line must lie in 
the region bounded approximately by 0,4 and 0 o 7 microns. 

(7) The absorbing substance must not be present in 

the atmosphere in sufficient quantity to appreciably af- 
fect the return signal. v 
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(8) The (half) half width of the line must be at 

o 

least of the order of 1 A in order that it include the 
entire spectral extent of the laser line. 

Liquids and solids, by virtue of their large inter- 
molecular forces, have absorption lines which are too broad 
for our purposes. Requirements (2) and (3) thus dictate 
the use of a gaseous absorber. The absorption lines of 
polyatomic molecules are always accompanied by a fine 
structure; * that is, by nearby secondary lines. Except in 
the case of H^, a molecule with extremely low mass, these 
nearby lines fall within 150 cm* 1 of the primary line, and 
thus render polyatomic gases undesirable due to require- 
ment (1). We thus limit our search to monatomic gases. 
Because of .the strong absorbing properties that we require, 
we limit our search to the resonance absorption lines of 
these gases. The only elements which produce stable mona~ 
tomic gases at reasonable temperatures, and whose resonance 
absorption lines lie within the wavelength region given in 
requirement (7) are the Group I (alkali) elements. Of 
these, the one that is to be preferred is sodium (The 
others are either too reactive, too rare, or both). 
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3. Characteristics of the resonance absorption lines of 
sodium* 

The Handbook of Chemistry and Physics indicates that 
sodium has two very strong resonance absorption lines at 
5890 A and 5896 A* Walsh (1961) gives Doppler (half) half 
widths of the 5890 A line as observed in flames at one 
atmosphere pi*cssaire and flame temperatures of 1000 K, 

2000°K and 3000°K. Penner (1968) shows that the Doppler 
(half) half width, is given by 

b D " C 1 T ^ (88) 

% 

and is independent of pressure. Figure 27 shows a curve 

of this form fitted to the data given by Walsh. Walsh 

also states that, in a sodium flame, where p « 1 atm and 

T ** 2200°K (broadening gas assumed to be N 0 ), the pressure 

broadening is approximately one half the Doppler broadening* 

From Figure 27, then, we determine that at 2200°K, 

b » 0.116 cm* 1 . Thus, b ** 0,058 cm" 1 at that texaper^ure. 
D P 

Pressure broadening, according to. Penner (1968), varies 
according to 

P . 

b 53 C 0 -r (89) 

p 2 

The pressure broadening (half) half width, b^, deduced from 


Doppler ( Holf) Half Width of Sodium Line of 5890 A , b (cm" 
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Walsh's data, is plotted as a function of pressure in 
Figure 28 for a temperature of 673 °K, the anticipated 

* • 

temperature of the sodium cell. 

/, ti - - .■ 

Fuchtbauer and Dinkelacker (1923) studied the shape of 
the mercury resonance line as a function ef the broadening 
gas and density. It is to be expected that, since. the 
resonance in mercury is not as intense as that in sodium, 
the mercury line is more readily pressure broadened, and 
thus ^ for mercury is greater than that for sodium at any 
given temperature and pressure. Thus, the results of 

it. . 

Fuchtbauer and Dinkelacker for mercury vapor, taken at 300 K, 
provide a means of checking the consistency^ of the results 

for sodium* Their figures, reduced to 673°K, are shown for 4 

• , 

comparison in Figure 28. The two sets of figures are 
Indeed compatible, with the values of for mercury great- 
er than those for sodium by a factor of about 2*0. For a * 

o 

sodium line width of \ A, these results Indicate that a ; . 
pressure of about 28 atm of N 2 is required. If argon is 
used as the broadening gas, the results of Fuchtbauer and 
Dinkelacker for mercury indicate that the results would be 
about the same. It should be noted herSTthat, at T * 673°K 
and p * 28 atm, .about 94 percent of the total line width 
is produced by pressure broadening. 




on 
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4. Resonance absorption cross section of the sodium mole - 
cule » 

The resonance absorption cross section of sodium was 
measured in the laboratory, using a Cary* Model 14 Spectro- 
photometer* A sample cell was constructed which gave an 
optical path length of about 2*5 cm (see Figure 29)* A 
small quantity of metallic sodium was introduced into the 
cell, and the cell was evacuated and sealed. By means of 
a heating coil wound around the cell, the sample tempera- 
ture was increased. The temperature of the cell was 
measured by means of an iron-cons tent an thermocouple with 
which it was in contact - 

The Handbook of Chemistry and Physics gives the vapor 
pressure of sodium corresponding to several temperatures 
above 712°K, If we assume that the equilibrium between 
metallic sodium and its vapor obeys the Claus ius-Clapeyron 
equation (Hess, 1959), that the vapor is an ideal gas and 
that the latent heat of vaporization varies linearly with 
temperature, then 

M(Na) db T 

+ in, — (90) 

R* dT T 0 

*Cary Instruments, 2724 South Peck Road, Monrovia, CA 91016* 


In 


e (Ha) 
e Q (Na) 



1 1 
T T. 
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* Sauereisen Insa-Lute Cement, No* 1* 


Figure 29. Construction details of temperature controlled 
sodium vapor absorption cell. 
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Using this equation and the data given in the Handbook, we 
find that 

h « 4.281 x 10^ erg-gm ^ 

dL . . J -I o -1 

•gj » 1.94 x 10 erg-gm - K 

Through use of these values, we may now obtain e(Na) for 
any temperature from (90). We obtain the number density of 
sodium molecules .in gaseous form at each temperature from 

n(Na) - SlM (91) 

XT 

Graphs of e(Na) .versus T and n(Na) versus T are given in 
Figures -30 and 31* 

The absorption of light passing through the sample cell 
° a 

at » 5890 A and at A 2 ~ 5896 A, as well as adjacent 

wavelengths, was measured for various cell temperatures 

through use of the spectrophotometer. Spectral slit width 
_ o 

was set at 0.5 A, corresponding to an anticipated laser 
bandwidth of that order of magnitude. Results of these 
measurements are shown in Figure 32. The transmission 
through the cell at the center of each line is given by 



Equilibrium Vapor Pressure Over Solid Sodium, e(No) (dyne- cm 


CM 



Figure 30. Equilibrium vapor pressure over solid sodium, e(Na) j 
c as computed using data given in the 
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Figure 31. Number density of sodium vapor in equilibrium with 
solid sodium# 



TRANSMISSION 
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5885 5890 5895 5900 5905 

* . 

WAVELENGTH (A) 


Figure 32. Observed transmission at various temperatures 
through absorption cell containing sodium vapor in equi- 
librium with solid sodium* 


0 * 014 * 
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- - r c 

= e 

(92) 

— e 

(93) 


The transmission at some wavelength, away from the 

influence of the two lines is given by 



(94) 


Through use of the third of these equations, we may elitni- ^ 
nate the unknown quantity / £' c from the other two, with the 
result; 


1 T* 3 

O' « In 

1 ”»« 'M 


(95) 




(96) 


Results cf the computations, shown in Table 11 and in 
Figures 33 and 34, indicate that the absorption cross sec 
tions vary somewhat with temperature (pressure). This is 
attributable to resonance effects as described by 
Fabelinskii (1968) for mercury vapor and other gases at 



Table 11 . Parameters of the absorption lines of sodium vapor at 5890 A and 5896 A, and factors used 
in their computation. 


T(°K) 

e(Na) from 
Eqh (90) 
(dyne*-OTi ) 

n fran 
Eqn (91) 
(csT*) 

T *1 

Tt 

A 3 

A 2 

T v * 
A 3 

(full) 

half 

width 

(full) 

half 

width 

AA 2 (2) 

Computed 
<5'x 1 (cm 2 ) 

Computed 
(cra ) 

1*7 9 

4.8 x 10“ 1 

7.3 x 10 12 

0.557 

0.590 

0.565 

0.590 

i .05 

0.975 

3.1 x IQ* 15 

2.5 X 10‘ ls 

9x7 

6.7x10° 

9.0 x 10 13 

0.503 

0.568 

0.523 

0.571 

1.13 

1 .09 

5.1* x 10 -16 

3.7 x 10‘ 16 

582 

2.8 x 10 1 

3.5 x 10 11 * 

0.1*15 

0.572 

0.1*58 

0.57U 

1.21* 

1 .21* 

3.7 x 10~ 16 

2.6 x 10" 16 

629 

1 .5 x 10 2 

1.7 x 10 15 

0.181 

0.517 

0.279 

0.533 

1.58 

1.65 

2.5 x 10' 16 

1 .6 x 10 

6?4 

5.2 x 10 2 

5.6 x 10 1 ; 5 

0.119 

0.436 

0.225 

0.1*50 

1 .28 

1 ,58 

9.3 x ID' 17 

5.0 x 10- 17 


* Transmission outside of the absorption lines appeared to have a slight tendency to increase with 
increasing wavelength (see ‘Figure 32). To ccmpensate for this, separate off-line transmission 
values were extracted for each line. 
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o Absorption cross section at 5890 A, 
x Absorption cross section at 5896 A t 


Temperature ( 


Figure 33. Computed absorption cross sections for the 
two resonance absorption lines of sodium vapor* 




Aitl uC.rtb!h , X r (Cm } 
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low pressures - 

For large optical thicknesses, we may neglect Z* c and 
write, for the center of the stronger absorption line: 



1 


^ T Al 


(97) 


Figure 34 shows the required path lengths for various 
combinations of center transmission and cell temperature, 
as obtained through use of this equation and the values of 


and n that were determined experimentally. The re- 


quired path length reaches a minimum near 650°K, since it 
depends on the quantity, 1/ which is a minimum at 


that temperature. 

Aside from the advantage of a relatively small cell 
si 2 e, we would prefer to operate our cell near this tempera- 
ture because in this region, cell absorption is relatively 
insensitive to small fluctuations in cell temperature * 
Suppose that we decide to operate the cell at 650°K, with 
a minimum transmission value of 10 5 . The required path 
length is 20 cm. 


5. Materials for construction of the cell* 

Having determined the size of the absorption cell and 
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the regions of temperature and pressure in which it is to 

operate, a search was conducted for materials that could 

possibly be used in its construction. 

Three samples of special alkali-resistant glass 

(Types 1720, 1723 and 1915) were obtained from the Corning 

o 

Glass Works* and were baked for about 70 hours at 650 K in 
a sealed chamber containing sodium. In all cases, a brown 
coloration, as well as severe etching of the surfaces was 
observed, rendering the glass totally unsatisfactory from 
an optical standpoint. 

The Westinghouse Electric Corporation**, a manufactur- • 
er of sodium vapor street lamps, was consulted. They 
stated that, while the glass used in street lamps is 
translucent, it cannot be made transparent. They suggested, 
however, that perhaps sapphire might fulfill our require- 
ments. Sapphire is a hard, inert , substance that is trans- 
parent to yellow light. 

Tyco Laboratories***, a manufacturer of industrial 
sapphire material was contacted. It was determined that 
they do not, at this time, manufacture optical quality 

* Mr. Robert Tichane, Corning Glass Works, Corning, NY 14830. 

** Mr. R.D. Hutchinson, Westinghouse Electric Lamp Division, 
•Bloomfield, NJ. 

***xy C o Saphikon Division, 16 Hickory Drive, Waltham, Mass 021 
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sapphire in pieces large enough for the desired purpose. 

6. Conclusions. 

It has been concluded that, while the use of a sodium 
absorption cell to selectively attenuate light at the 
transmitted frequency has merit, it is not feasible to 
implement such a system at this time, due to the tempera- 
tures and pressures involved, as well as the high chemical, 
activity of sodium. Perhaps with the development of 
tunable lasers which operate efficiently with a narrower 
spectral output, it will not be necessary to broaden the 
sodium line as much, thereby increasing the practicability 
of a system of this sort. 

Use .of an absorption type filter is particularly ap- 
pealing in that it all but does away with the stringent 
optical collimation requirements inherent in all inter- 
ferometric devices. 
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. Chapter V. 

SIGNAL TO NOISE RATIOS 

List of Symbols 

Effective collection area of receiver (cm 2 ). 

Speed of light (cra-sec* 1 ) • 

Fraction of available Raman energy used by filter 
with (half) half width, h. 

, "*1 \ 

Filter (half) half width (cm- ). 

Relative, background current received in Raman . 
channel. 

Relative signal current received in Raman channel. 
System constant. 

Number of laser pulses per observation. 

Number density* of air (N 2 and 0^) molecules (cm- )« 
Electronic charge (coul). 

, *»1 , 

Responsivity of photomultiplier tube (tnA—W ). 

r -2 , -lv-i _ ~n . 

Background radiance L W-cra -(cm ) -ster - J . 

Transmission of filter at center of its transmission 
curve. 

transmission at frequency ^ of filter centered 

at V * with (half) half width, h. 
c 

Atmospheric transmission -* molecular component. 
Atmospheric transmission - particulate component. 
Transmitted energy (J). 

Height (cm). . 



Bandwidth of signal processing apparatus (sec ). 

Wavelength ) . 

, .-l x 

Frequency (cm }• 

Frequency at^sjrti^^ filter transmission curve is 
centered (cm ). 

Signal to noise ratio. 

Effective Raman backscatter cross section of air 
(cni -ster** 1 ) . 

Receiver acceptance angle (sfcer). 
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A, Absolute signal to noise ratios at 0,1 M . 

In Equations (2) through (4), it was shown that the 
signal current received in the Raman channel of a mono- 
static pulsed lidar system may be written as 

i s (h,z) = k r e'j.nC jf m QT p ) 2 ( 98 > 


The system constant k r is given by 



W t A r RcE(h)T c 



(99) 


In the following discussion, T c will be assumed to be in- 
dependent of the type of filter employed. This Is approxi- 
mately true for relatively broad filters. For very narrow 
filters,. the attainable center transmission can be expected 
to decrease somewhat with decreasing width. 

The current received due to background radiance is 




i b < h > 


9$RA r m c 1 T y Ch,V c )d()/-V c ) (100) 


y~ca 


For a Gaussian filter function. 



(101) 
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and 


r" ' 

/v-V\ 2 


\ ' 

-(in 2)1— g—J 


\T*<h,v c ) d <v -y) - 

i f * V 

e V h ' - h 

In 2 



wm. 


( 102 ) 


Thus* 


i b ( h ) 


^RA r UT c h 



In 2j 


( 103 ) 


If both the signal and background levels were of constant 
Intensity* the observation would be noise-free* In that the 
constant background could be precisely determined by monitor- 
ing the current received with the laser not in operation* and 
subtracted out during the computations. In reality, both the 
signal and background levels are subject to random statistical 
fluctuations about their mean values* even though the means 
themselves may be constant « If we treat the arrivals of 
photons as rare events, then we may apply Poisson statistics 
and assume that the RMS fluctuation in the number of photons 
received by a receiver with bandwidth ^f is equal to the 
square root of the mean number received* It is this random 
fluctuation about the mean which constitutes Che noise. 

Thus* the signal to noise ratio is given by 


£ (h* z) 


i s (h,z) 


~^2q At i g ( h , 2 ) + i^Chj 


( 104 ) 
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We can evaluate (h,z) approximately by considering a 
model atmosphere and a typical laser radar system. We 
choose number densities given in the U» 5 . Standard Atmos* * 


>here» 1962, and the following system parameters; 



In the above discussion we have assumed that, in collecting 
the signal and background photons, the laser is fired only 
once* Signal to noise ratios may be improved somewhat by 
averaging over a number of laser pulses, if the pulses are 
sequenced in such a way that 

1) consecutive pulses are sufficiently far apart that 
the fluctuations about the means are independent for each 
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pulse, and 

2) the total time period involved in taking the 
observation is short enough for the mean value to be considered 
as stationary. 

When the above requirements are met, the signal to noise 
ratio increases as Clearly, the times involved in 1) and 

2) above depend upon the temporal, spatial and spectral scales 
of the observation 

For the sake of simplicity, and due to the fact that 
one hundred pulses are required to produce a factor of ten 
increase in the signal to noise ratio, we shall pursue our 

♦ 

analysis with the assumption of a single pulse observation. 

Results of the computations, shown in Figures 35 and 36, 
indicate that, for nighttime work, the signal to noise ratio 
is shot-noise limited; that is, 

« S 

and thus, the signal to noise ratios are proportional to 
1/E(h). Since 

lim E(h) = 1 

this dependence results in a slow increase in £ at 
moderate half widths, reaching an asymptotic value for large 
widths . 

For daytime work, i s and are of the same order of 
magnitude at low altitudes, with i fe predominating over i g 
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Figure 35. Nighttime signal to noise ratios ^ 

radar system, showing the effect of filter half wine 
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at greater heights* In the background-limited case, when 

* b » l . 

the signal to noise ratios are proportional to E(h)/Vh * 
Here, if E(h) Increases more rapidly with h than doesYh" , . 
the signal to noise ratio increases with increasing h. 

If E(h) increases less rapidly with h than doesVTT , then 
the signal to noise ratio decreases with increasing h. For 
purposes of analysis, we wish to locate the point which 
gives the optimum signal to noise ratio for a background- 
limited signal; that is, the point where 


or 


d E(h) 
dh ViT 



d E(h) 

— E(h) » % 

dh h 


(105) 


(106) 


For both the 0- and S- branches, this point lies near 

h ** 50 cm \ For signal to noise ratios falling between 

the background-limited and shot-noise limited values, the 

optimum half widths are shifted toward larger values. 

For the system parameters assumed here, the optimum 

—1 

half width falls at soma value greater than 70 cm , but 

the signal to noise ratio does not deteriorate significant- 

~1 

ly as long as h I 3 greater than about 30 cm . 
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B, Determination- of wave len gth, of; optimum signal to noise 
ratio, 

A computation was .aad^* of signal to noise ratio as a 
function of wavelength (daytime operation) , with a view 
toward finding the opting wavelength regicr: in which to 
operate, and tc determine the penalty paid for operating 
away from thtt wavelength. 

Table 12, shows the magnitude of each of the wavelength- . 
dependent variables affecting the signal levels, relative 
to their magnitudes at 0.7a* * in the verge 0*3// ^ 0«?&, 

end Table 13 *hrv?s the variation of background level, i b j wit 
wavelength. The signal level, i g , is also a function of 
i/ss , so the optimum wavelength may also be expected to be 
a function -of height, 

Relative signal to noise ratios derived using these 
data s together with the complications of the previous chap- 
ters, are showv in Figure 37, The ratios are shown to be 
optimally large near about 0.^5// 3 but the ratios are 
shown to wry by i. c j * than a factor of three between 0,35,4/ 

and 0.?// 
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Figure 37. Relative signal to noise ratios for various heigh 
showing the effect of varying the laser wavelength. 
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Table 12* Magnitude (relative to that at 0.7,#) of each of 

the wavelength-dependent variables affecting the signal 

level, i . 

* s 


AU) 

t 

1 

A 4 

O ' 2 

After 

Elterman 

(1968) 

R 

After RCA 
Corp. (1970); 

j 

w t 

(product of 
preceding 
four 
columns) 

0.3 

29. 

7.6 x 10“ 5 

1.09 

0.10 

2.5 x 10~ 4 

0.4 

9.4 

0.43 

1.78 

0.25 

1.8 

0*5 

3.8 

0.72 

1.55 

0.50 

2.1 

0.6 

1.9 

0.83 

1.36 

0.75 

1.6 

0.7 

1.0 

1.0 

1.00 

1.0 

1.0 


Table 13* Magnitude (relative to that at 0.7^) of each of 
the wavelength-dependent variables affecting the back- 
ground level, i^. 


A (A) 

n 

R 

After RCA 
Corp. (1970) 

i 6 

. (Product of preceding 
two columns) 

0.3 

0.14 

1.09 

0.15 

0.4 

0.28 * 

1.78 

0.50 

0.5 

0.46 

1.55 

0.71 

0.6 

0.69 

1.36 

0.9V 

0.7 

1.0 

1.00 

1.0 


.i* 




166 


Chapter VI. 

CONCLUSIONS 

The conclusions to be drawn from this research may be 
divided into two parts: those pertaining to the theoreti- 

cal aspects of the technique and those pertaining to its 
practical application. 

It has been shown that the rotational Raman spectrum 
of light scattered from air contains regions in which the 
effects of both temperature and particulate content are 
negligible. This fact is of importance not only for the 
present study,' where it is desired to extract turbidity 
data from the elastic scattering, but may also be helpful 
to investigators such as Salzman, et.al. (1971), who are 
attempting to utilize the temperature-sensitive portions 
of the rotational Raman spectrum of air to obtain remote 
temperature soundings. 

It has been established here that particulate effects 

o 

are negligible at 4380 A, the operational wavelength of 
the Raman spectrophotometer. It is quite possible, how- 
ever, that off-frequency scattering from particulates could 
become significantly large at shorter wavelengths, where 
resonance scattering and flourescence become important. It 
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is well known that resonance scattering and fluorescence 

cross sections may be orders of magnitude greater than 

those discussed here. In any case, it was shown "in this 

o 

study that wavelengths below about 4500 A are not desire 
able for turbidity measurements because of the large rate 
of change of atmospheric transmission with wavelength. 
However, this caveat should be borne in mind by any in- 
vestigator contemplating the extension of these results 
into the blue or ultraviolet portion of the spectrum. 

With the attainability in recent years of high resolu- 
tion, low noise rotational Raman spectra, the rotational 
Raman spectrum has been routinely observed and, indeed, is 
even used in the alignment of Raman spectrophotometers. 

By utilizing some of these observations, an attempt was made 
in this study to check for consistency the polarizability 
anisotropies for and 0^ obtained by earlier investigators. 
While there is still some doubt as to the absolute values of 
the anisotropies, it can at least be said that the values for 
and 0 9 obtained by Bridge and Buckingham (1966) are 
correct relative to each other. 

Investigation of the practical applicability of this re- 
mote sensing system reduced primarily to the development of 
a technique for spectrally and spatially separating the 
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Raman scattering from the much stronger elastic scattering. 
This phase of the study gave the author considerable in*" 
sight into the problems involved in optical signal proces- 
sing. Most importantly, the computations revealed (to the 
author, at least) the basic unity of all interferometric 
systems- the problem of focussing rays of light through a 
slit on to the face of a diffraction grating is transformed 
to a problem of similar difficulty if one decides instead 
to collimate the rays and use an interference filter. 

It was ascertained that the rigid requirements for an 
interferometric system could be relaxed somewhat if an 
absorption cell could be used to preattenuate the elastical- 
ly-scattered light. Such a cell, if it were to be effective, 
however, would present construction problems distinct from, 
but perhaps more difficult than those encountered in an 
interferometric system. 

All told, it appears that, although a suitable filter 
system is within the . realm of possibility today, the present 
state of- the art makes such a system expensive and/or inef- 
ficientV Despite these inefficiencies, it is found that a 
typical laser radar system can produce daytime signal to 
noise ratios of the order of 10 db to heights of at least 



5 knu For nighttime work, 10 db signal to noise ratios 
are achievable to heights of at least 15 km. 
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Appendix 

1 0 COMPUTATION OF THE RATIO y 2 (Nj/ y 2 (0 ~ ) FRQM 0BSERVE ^ 
DATA* . 

The intensities of Individual nitrogen and oxygen rota- 
tional Raman lines, relative to a common datum level, are 
given by 

Iy (N 2 ) 
ly (0 2 ) 



wv * z 


(107) 


W°2> + 2 


(108) 


where Z is the zero displacement of the scale. If we recog- 

2 

nize that intensities are proportional to nJ , • then we 
may write 

Kn(N 2 )S 2 (N 2 ) + Z (109) 

* . * 

Kn(0 ? )if 2 (0,) + Z (110) 

Thus, if we plot a graph of IyCNj), the intensities 
extracted from the results of Miller, et.al. (1968), versus 
[l y (N 2 ) i I xot^ N 2' , -l J the resulting curve should be a straight 


ly(0 2 ) = 


I x y (i*«) 


^TOt^V 


X y (N,) = 




<N 0 ) 


WV 
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line, with slope proportional to 

Kn(N 2 ) JT 2 (N 2 ) 

The analogous graph for has a slope proportional to 

Kn<0 2 )* 2 (0 2 ) 

The constant of proportionality, K, is of course, the 
same for both curves, as is the intercept point with 
the coordinate axis (the 2 ero displacement of the scale). 

The two curves plotted in this manner are shown in Figure 
38, together with the best fit straight lines through the , 
points , 

It should be- noted that ly /l Xot is a function of 
temperature. According to Miller (1970), their experiment 
was conducted at, or slightly above room temperature, but 
no attempt was made to measure the temperature. Graphs 
were plotted for several temperatures, and the best straight 
line fits were obtained for T = 287°K C 

The quotient of the two slopes was found to be 

n(N,)* 2 <N 7 ) 

= ~ = 1.496 

n(0 )iT(0 ) 



0 I 2 3 4 5 6 7 8 9 10 It 

[I„/I tot 3 ( Calculated ) 

Figure 38, Intensities of lines in the rotational Raman spectra of and 0 2 in air, h-j 
as observed by Miller, eto‘al 0 (1968), versus intensities computed in Chapter III* 
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2 * TESTING OF THE INTERFERENCE FILTER* 

It is desired to determine whether the interference 
filter for the Raman channel has a transmission of 
10” 4 to 10~ 5 as specified by the manufacturer- Fur- 
thermore , it is desired to examine the sensitivity of 
this transmission value to the angle of the incident 
light. 

An optical system was constructed to determine the 

o 

filter parameters at 6943 A for a range of incidence 
angles near zero degrees « A schematic diagram of the 
device is shown in Figure 39 . The mirror M, actually an 
unsilvered piece of glass, has a low reflectance, and is 
used to direct a small percentage of the laser output 
into the test apparatus. The lens combination and 
was chosen to produces 

lo A collimated beam of light with a diameter about 
equal to that of the interference filter (one inch), 

2. A sensitivity of the degree of collimation to the . 

location of That is, it was desired to vary the beam 

•2 

width of the light striking the filter from about one 
degree of convergence through about one degree, of divergence 

3, A spot of light nearly constant in size at the 
location of F, the interference filter, regardless of the 





Figure 39... Schematic diagram of device constructed to 
test transmission characteristics of interference 
filter. •./ ; \ • 
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location of L^. 

is used to concentrate the light transmitted through 
the interference fitter on to a cadmium sulfide photocell, P. 
The iris diaphragm, X, is used to adjust the size of the 
spot* 

The lens L 2 is located on a vernier mount which allows its 
position along the optical axis to be adjusted over about 
±2*5 cm from its mean position. The interference filter, F, 
is mounted’ in such a way that it can be tilted up and down 
and from side to side* 

The entire apparatus is enclosed in a lightproof box, 
which is painted on the inside with flat black paint to mini- 
mize diffuse reflection* 

The device is designed so that the interference filter, 

F, can readily be removed and replaced by broad band absorp- 
tion type filters whose transmissions are known, and which 
lie in the same range as that of the interference filter. 

Calibration of the system and testing for photocell 
linearity was accomplished by measuring the photocell out- 
put with filter combinations as shown in Table 14 and in 
Figure 40. 

Calibration of the vernier moxint supporting the filter 
consisted of replacing the filter, F, with a plane mirror. 
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Figure 40, Results of test for linearity of photocell output 
when operating in the configuration shown irj Figure 39, 



and applying trigonometry to the location of the reflected 
light to obtain the mirror angle. Results of these measure- 
ments are shown in Figures 41 and 42. 

Figure 43 shows resulting transmitted intensities obtained 
when the interference filter was mounted at various vernier 
settings. Transmitted intensity through a neutral density 
filter with a transmittance of 2 x 10~ 5 is shown for compari- 
son. Several pulses of the laser are recorded for each 
setting in order that a suitable average reading could be 
obtained. The variability from pulse to pulse is attributable 
to the laser, the output of which fluctuates considerably. 

Figure 44 shows the average intensity values plotted 
against the vernier (angular) setting. When the filter is 
optimally tilted (at some angle very close to zero degrees) , 
transmittance s as low as about 5 x 10 ^ are attainable. 

For off-angle radiation, the transmittance rapidly increases 
reaching nearly 2 x 10 for incidence angles cf less than 
one degree. 
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Figure 44. Intensity values of Figure 43. plotted against incid 
angle deduced from vernier setting. Relative transmittance of 
the neutral density filters is shown for comparison. 



